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EXECUTIVE  SUMMARY 


The  Free  Electron  Laser  (FEL)  weapon  is  excellent  for  providing  state-of-the-art 
air  defense.  It  is  also  effective,  if  need  be,  in  offensive  operations.  It  is  possible  for  this 
weapon  system  to  achieve  a  soft  kill  whenever  needed  by  adjusting  lethality.  When 
combined  with  its  accuracy  makes  it  a  valuable  tool  in  operations  where  the  requirement 
is  not  to  only  win,  but  also  to  avoid  casualties. 

One  transitional  step  toward  the  1  MW  output  power  required  for  a  laser  weapon 
is  the  proposed  100  kW  upgrade  of  the  Thomas  Jefferson  National  Accelerator  Facility’s 
FEL.  At  present,  the  TJNAF  FEL  is  the  most  powerful  FEL  in  the  world  and  operates  at 
more  than  2  kW  average  power. 

Our  multimode  simulations  showed  that  this  development  is  feasible,  and  for  the 
current  TJNAF  lOOkW  FEL  parameters  with  <2  =  4.2,  the  inversely  tapered  undulator 
A  =  -7t  achieved  the  highest  final  power  at  a  desynchronism  of  d  =  0.04 .  The  output 
power  produced  was  P  -120  kW,  which  is  above  the  lOOkW  objective.  The  small  values 
of  the  induced  electron  energy  spread  (-3.5%)  made  it  possible  to  explore  a  design  with 
lower  10%  losses  of  (Q  =10)  in  an  optical  cavity.  The  higher  Q  factor  was  chosen  to 
increase  the  efficiency  of  the  system,  while  maintaining  the  induced  energy  spread  well 
below  the  15%  engineering  upper  limit.  Simulations  results  with  Q  =10  design  indicate 

that  the  inverse  step- taper  undulator  A  =  —%  achieved  the  highest  final  power 
P  =  190kW  at  a  desynchronism  of  d  =0.01.  The  inverse  linear  taper  8  =  -271  achieved 
the  highest  final  power  P  =200kW  at  a  desynchronism  of  2=0.01.  For  all  the 
aforementioned  cases  the  steady- state  gain,  the  final  steady- state  power,  and  the  induced 
electron  spread  as  a  function  of  desynchronism  and  taper  rates  were  determined. 

The  validity  of  our  results  was  verified  against  experiments  conducted  in  the 
TJNAF  FEL  facility  at  lower  power.  FEL  operation  as  a  function  of  various  taper  rates 
and  desynchronism  values  was  studied.  The  purpose  of  the  simulations  was  to  compare 
experimental  observations  with  theoretical  analysis,  and,  by  using  a  wider  range  of 

xvii 


parameters  than  allowed  in  the  experiment,  to  extract  more  physical  meaning  from  the 
results. 


The  simulations  and  the  experimental  data  agreed  well  and  for  the  most  part  were 
consistent  with  analytic  theory.  This  theory  indicates  that  as  the  taper  increases,  the 
desynchronism  width  decreases.  The  experiment  shows  this  trend,  but  not  as  clearly  as 
the  simulation  data.  Only  one  point  shows  a  significant  deviation  between  the  simulation 
and  experiment.  This  has  resulted  in  more  examination  of  the  experimental  data, 
revealing  that  the  conditions  changed  during  the  measurement  of  the  point  in  question. 

The  maximum  gain  with  no  tapering,  5=0,  occurred,  as  would  be  expected.  As 

the  taper  rate  5  increases  or  decreases,  the  gain  decreases  from  160%  to  about  40% 
symmetrically  around  5=0.  At  positive  or  negative  taper  rates  near  5  =  ±871  ,  the  gain 
plateaus  observed  decreased  much  more  slowly  with  an  increasing  taper  rate.  As  the  taper 
rate  increases  in  magnitude  from  |5 1  =  8tc  out  to  |5 1  =  2471  ,  the  gain  decreases  from  about 
40%  down  to  about  10%.  The  gain  plateau  at  large  values  of  taper  begins  at  a  value  of 
taper  5  that  causes  the  gain  spectrum  to  change  shape  so  that  there  are  two  peaks  with 
nearly  the  same  peak  gain.  As  the  taper  is  increased  and  the  peak  available  gain 
decreases,  the  gain  spectrum  acquires  multiple  peaks  with  comparable  gain.  As  the 
primary  gain  peak  decreases,  other  surrounding  gain  peaks  increase  in  comparison.  This 
can  also  occur  in  the  untapered  case,  but  results  in  a  smaller  gain  peak.  As  the  peak 
available  gain  increases,  or  decreases,  at  any  value  of  tapering,  the  desynchronism  curve 
width  A d  increases  or  decreases  correspondingly.  With  more  gain  available  there  are 
more  values  of  desynchronism  that  are  above  the  threshold  which  makes  the 
desynchronism  curve  wider. 
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I.  INTRODUCTION 


Directed  energy  weapons  (DEW)  is  not  a  new  initiative  that  “just  appeared.”  It  is 
the  consequence  of  changes  in  political  philosophy  and  technology.  The  emergence  of 
more  and  more  power  sources  of  microwave  energy  is  clearly  a  key  factor  in  their 
applicability.  Furthermore,  many  of  these  sources  are  available  at  reasonable  cost  and 
satisfactory  sizes. 

DEWs  have  been  on  the  on  the  U.S.  military’s  “love  to  have”  list  for  several 
decades.  Within  the  United  States,  all  three  services  have  been  actively  developing  laser 
technology.  One  of  the  most  impressive  land-based  Tactical  High-Energy  Lasers  (THEL) 
is  the  “Nautilus”.  It  is  a  joint  U.S.  Army  and  Israeli  project,  which  uses  a  Mid- Infra-Red 
Advanced  Chemical  Laser  (MIRACL).  At  present,  it  is  being  tested  at  White  Sands 
Missile  Range  in  New  Mexico  and  has  been  successfully  tested  against  the  supersonic 
BM-21  122mm  artillery  rocket.  The  Nautilus  program  is  still  being  tested  and  reached  a 
milestone  when,  in  June  2000,  a  test  firing  resulted  in  the  destruction  of  a  Katyusha  type 
missile  in  flight  [Ref.  1], 

The  U.S.  Air  Force,  as  well  as  the  U.S.  Army,  in  conjunction  with  the  Air  Borne 
Laser  (ABL)  program,  adapted  the  chemical- oxygen- iodine  laser  (COIL)  to  a  flyable 
model  against  missiles.  This  project,  which  was  on  schedule  and  within  budget,  is  now 
transitioning  from  the  early  concept  design  to  operational  use. 

The  U.S.  Navy  is  also  interested,  and  on  the  24th  of  April  2001  the  Commander  in 
Chief  of  the  U.S.  Atlantic  Fleet  sent  the  following  letter  to  the  Chief  of  Naval  Operations: 

...and  has  proven  that  Speed  of  Light  weapons  can  be  very  effective 
against  these  small  high  speed  threats.  Such  a  laser  weapon  would  offer 
our  Naval  forces  an  extremely  versatile  weapon  to  counter  numerous  soft 
and  hard  targets.  A  High  Energy  Laser  weapon  can  be  designed  to  deliver 
energy  that  can  track,  warn,  damage,  mission  kill,  and  if  need  be,  destroy  a 
threat.  I  believe  it  is  exactly  this  type  of  weapon  system  that  our  forces 
need  in  the  littoral  environment  where,  even  though  the  threat  may  not 
always  be  as  sophisticated  as  a  highly  maneuverable  cruise  missile, 
intentions  are  often  more  difficult  to  discern  and  timelines  are  extremely 
short . 
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A  Free  Electron  Laser  (FEL)  as  a  weapon  is  an  excellent  approach  to  fulfilling 
present  and  future  requirements.  It  could  provide  state-of-the-art  air  defense  not  only  to 
ships  but  also  to  sensitive  areas  such  as  Naval  and  Air  Force  Bases.  It  is  capable  of 
achieving  a  high- probability  hard  kill  against  all  cruise  missiles.  It  increases  the 
engagement  range  and  minimizes  engagement  time  and  thus  reduces  the  threat  of 
simultaneously  arriving  cmise  missiles.  Reliability  is  almost  100%.  Since  only  light  and 
vacuum  are  involved  in  the  weapon  system,  minimal  maintenance  is  required.  Its  infinite 
magazine,  its  instantaneous  reaction  at  the  speed  of  light  time  and  the  rapid  re¬ 
engagement  point  offer  many  advantages  over  gun  and  missile-based  close-in  weapons. 

It  is  also  possible  for  this  weapon  system  to  achieve  a  soft  kill  whenever  needed 
by  adjusting  the  lethality.  The  latter,  combined  with  its  pinpoint  accuracy,  makes  it  a 
valuable  tool  in  operations  where  the  requirement  is  not  to  only  win,  but  also  to  avoid 
casualties  [Ref.  2].  This  fact  is  present  in  all  peacekeeping  operations,  as  it  was  in 
Kosovo  and  now  in  Afghanistan.  This  philosophy  has  been  another  push  toward  the 
deployment  of  DEWs  and  FELs. 

Presently,  FEL  research  is  widespread  and  advancing  on  many  fronts.  The 
impetus  comes  from  some  unique  advantages  the  device  offers  over  other  types  of  lasers. 
[Ref.  3]: 

•  “Potential  for  producing  very  high  power  by  extending  technology 
developed  for  existing  electron  accelerators.  The  evolution  of  power 
densities  in  magnetrons,  klystrons  and  gyrotrons  have  reached  their  limits 
while  the  limit  for  FELs  have  not  yet  been  reached  as  Figure  1  indicates. 
There  the  product  of  the  average  power  and  the  square  of  the  radiation 
frequency,  as  a  figure  of  merit,  is  shown.  These  available  microwave 
powers  were  unthinkable  just  a  few  years  ago. 

•  Prospects  for  wall-plug  efficiency  more  than  10  percent  or  more  at  high 
powers. 

•  Potential  for  tunable  operation  from  the  millimeter-wave  region  to  the 
extreme- ultraviolet  or  soft  x-ray  region,  although  no  single  device  would 
operate  over  such  a  broad  range.  (The  broad  tuning  range  comes  from  the 
fact  that  the  electrons  are  not  bound  when  they  emit  light;  it  would  be 
achieved  by  varying  the  electron  energy  and  the  spacing  and  strength  of 
the  magnetic  field)" 
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Figure  1.  Progress  in  the  Development  of  a  Variety  of  Vacuum  Electronic  Sources  of 

Coherent  Radiation  as  Measured  by  the  Evolution  of  the  Product  of  the  Average  Power 
and  the  Square  of  the  Frequency.  “From  [Ref.  4].” 

Chapter  II  of  this  thesis  describes  the  FEE  weapon  concept  and  the  improvements 
gained  in  the  field.  Present  limitations  mentioned  at  the  Workshop  2001  in  Virginia, 
USA,  compare  the  FEE  with  other  weapons  (HEL  and  conventional)  and  address  its 
benefits. 

Chapter  III  gives  an  overview  of  the  theoretical  background  and  discusses  the 
physics  behind  the  FEL. 

Chapter  IV  describes  the  results  of  simulations  for  the  proposed  100  kW  upgrade 
to  the  TJNAF  FEL.  This  upgrade  is  a  transitional  step  for  the  1  MW  output  power.  This 
work  was  originally  presented  at  the  23rd  International  FEL  Conference  in  Darmstadt, 
Germany  and  has  been  accepted  for  publication  in  Nuclear  Instruments  and  Methods  in 
Physics  Research. 

Chapter  V  presents  the  analysis  of  simulation  results  based  on  experiments 
conducted  at  TJNAF.  Simulations  and  experimental  results  are  compared  to  analytical 
theory. 
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n.  FEL 


The  FEL  is  a  device  which  extracts  light  energy  from  a  beam  of  relativistic 
electrons  passing  through  a  spatially  periodic  magnetic  field.  John  Madey  first  introduced 
the  concept  in  1970  [Ref.  5].  It  is  unlike  other  lasers,  such  as  solid-state,  semiconductors 
and  liquid  lasers,  because  it  relies  cn  light  emitted  by  electrons  that  are  not  bound  to 
atoms. 

An  FEL  consists  of  an  electron  beam  source,  a  periodic  transverse  magnetic  field 
(an  “undulator”  or  “wiggler”  magnet),  and  an  optical  resonator  as  shown  in  Figure  2. 

Electron  beam 


resonator 


mirror 


Figure  2.  The  FEL  Schematic. 


The  undulator  magnet  in  an  FEL  imposes  a  transverse  acceleration  on  the 
electrons,  resulting  in  the  generation  and  amplification  of  light,  which  co-propagates  with 
the  electron  beam  through  the  magnet.  As  in  a  conventional  laser,  the  resonator  mirrors 
provide  feedback  around  the  amplifying  medium,  permitting  the  stored  radiation  to  build 
up  to  saturation  and  produce  a  very  powerful,  coherent  laser  beam. 

A.  WEAPON  CONCEPT -PRINCIPLES  INVOLVED 
1.  The  Way  the  Target  is  Destroyed 

A  laser  beam  in  the  basic  fundamental  mode  closely  approximates  a  Gaussian. 
The  intensity  drops  off  from  a  maximum  and  has  no  side  lobes.  The  typical  laser  beam  is 
quite  narrow  at  no  more  than  a  few  millimeters  in  diameter.  Furthermore,  when  operating 
at  that  mode,  the  minimum  waist  of  the  optical  beam  is  where  the  highest  energy  density 
occurs. 
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Since  the  beam  resembles  a  truncated  plane  wave,  it  is,  of  course,  highly  spatially 
coherent.  Its  directionality  may  be  thought  of  as  a  manifestation  of  that  coherence.  Laser 
light  is  quasimonochromatic,  and  generally  has  an  exceedingly  narrow  frequency 
bandwidth  indicating  that  it  is  highly  temporally  coherent.  Another  attribute  is  the 
amount  of  radiant  power  that  can  be  delivered  in  that  narrow  frequency  band.  The  laser  is 
distinctive  in  that  it  emits  all  its  energy  in  the  form  of  a  narrow  beam.  In  contrast,  a  100W 
incandescent  light  bulb  may  give  out  considerably  more  radiant  energy  than  a  lower- 
power  cw  laser,  but  the  emission  is  incoherent,  spread  over  a  large  solid  angle  with  a 
broad  bandwidth. 

A  good  lens  can  totally  intercept  a  laser  beam  and  focus  all  of  its  energy  into  a 
minute  spot  whose  diameter  varies  directly  with  wavelength  X  .  Spot  diameters  of  just  a 
few  thousandths  of  an  inch  can  readily  be  attained  with  lenses  that  have  a  conveniently 
short  focal  length.  The  laser  intensity  that  can  readily  be  generated  in  a  focused  beam  is 
over  1017  W/cm2,  in  contrast  to  an  oxyacetylene  flame  having  an  intensity  of  roughly  103 
W/cm".  A  focused  laser  beam  of  a  few  kilowatts  can  bum  a  hole  through  a  quarter- inch 
stainless  steel  plate  in  about  10  seconds.  [Ref.  6] 

Studies  and  scaling  experiments  have  shown  that  a  liter  of  a  missile’s  material  can 
be  destroyed  in  a  few  seconds  with  1  MW  of  laser  power  [Ref.  7].  One  liter  corresponds 
to  a  hole  with  the  dimension  10cm  x  10cm  x  10cm  which  is  generally  enough  to 
structurally  disable  the  missile.  Missile  destruction  can  be  achieved  in  different  by  an 
explosion  of  the  missile  warhead,  or  an  aerodynamic  instability  causing  the  missile  to 
breakup. 

A  targeting  system  on  a  launching  platform  will  control  the  beam  directed  to  the 
missile  and  bum  a  hole  in  it.  Unlike  conventional  bullets  or  missiles,  light  instantly 
travels  to  the  target  and  does  not  suffer  from  gravitational  effects.  The  Acquisition 
Tracking  and  Pointing  (ATP)  system  is  in  use  and  has  been  successfully  tested  in  the 
MIRACL  project  many  years  ago  [Ref.  1].  Acquisition  includes  the  detection  of  the 
target  in  the  tracking  system  using  infrared  radiation.  Tracking  begins  when  a  series  of 
consecutive  observed  positions  collected  to  allow  for  the  filing  of  the  azimuth,  elevation 
and  time.  Connecting  these  points  permits  the  construction  of  nominal  track  data  that 
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point  the  telescope  in  order  to  keep  the  image  of  the  target  centered  in  the  tracking  focal 
plane.  The  system  temporal  response,  which  includes  a  camera  processor  and  motors, 
must  be  fast  enough  so  that  the  target  does  not  disappear  from  the  camera  field  of  view. 
The  pointer/tracker  will  be  stabilized  against  unintended  motion  by  using  a  high  speed 
gyroscope  and  inertial  measurement  system.  ATP  could  cost  $6  million,  including  the 
multiple  target  acquisition  sensor  and  stabilization  system  [Ref.  8]. 

2.  The  Threat  -  Requirements  from  a  Weapon 

Targets  are  evaluated  with  respect  to  their  maneuverability,  size,  and  vulnerability 
to  attack.  Target  maneuverability  imposes  a  requirement  that  the  missile  must  be  capable 
of  terminal  maneuverability  to  effect  lethality. 

The  airborne -target  spectrum  is  quite  varied  and  consists  of  low  altitude,  slow- 
moving  targets  in  ground  clutter,  high- altitude  fast- moving  targets,  head-on  targets  with 
tremendous  closing  speeds  and  side-encounter  targets  with  high-g  turning  requirements. 
If  the  threat  is  other  than  airborne-targets,  no  known  limitations  are  implied.  In  contrast, 
the  surgical  accuracy  of  the  system,  the  absence  of  fragments  and  the  instant  battle 
assessment  make  the  system  operable  for  many  targets  and  environments.  Range  and 
time-to-range  is  a  factor  that  does  not  greatly  affect  the  FEL  weapon  design.  The 
engineer  need  not  consider  design  features  affecting  range  and  time-to-range  such  as 
drag,  wings,  controls,  lifting  surfaces,  and  other  parasitic  losses.  The  incredibly  small 
time  of  flight  of  33  (isec,  or  the  time  it  takes  the  beam  to  reach  the  missile  at  10  km,  is  the 
key  issue.  However,  environmental  constraints  and  atmospheric  transmittance  have  to  be 
taken  into  account. 

Cost  per  shot  must  be  much  less  than  cost  of  threat  target,  with  high  lethality  per 
shot.  For  example,  an  Unmanned  Aerial  Vehicle  (UAV),  which  is  used  for 
reconnaissance,  may  be  worth  less  than  an  expensive  sea-sparrow  missile.  The 
motivation  for  cheap,  effective  defense  does  not  imply  that  we  are  not  going  to  use  all  the 
means  of  defense,  whenever  human  lives  are  placed  in  danger.  Additionally,  the 
individual  engagement  time  must  be  short  in  order  to  engage  a  large  number  of  targets.  A 
dwell  time  of  3  seconds  makes  it  possible  to  handle  several  incoming  supersonic  targets 
simultaneously. 
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3.  FEL  Comparison  with  Existing  Weapons 

Present  systems  detection  is  limited  to  the  horizon  for  low  level  flying  threats. 
Unlike  the  antiballistic- missile  kinetic  kill  systems  designed  for  missile- terminal  phase- 
defense,  the  FEL  mode  of  attack  offers  several  key  advantages.  First,  the  target  is  much 
slower  moving  and  there  is  less  urgency  for  targeteers  to  worry  about  distinguishing 
decoys.  Additionally,  a  significant  advantage  of  the  FEL  is  that  it  is  environmentally 
friendly  as  it  does  not  generate  the  toxic  effluents  common  to  other  types  of  weapons. 
This  issue  of  “the  exhaust  gas  contaminants”  is  serious  since  the  crew  could  be  harmed, 
and  equipment  damaged.  Furthermore,  another  disadvantage  is  the  need  for  a  logistics 
supply  of  hazardous  and  toxic  chemicals.  The  latter  probably  is  the  strongest 
disadvantage  since  no  one  today  wants  to  handle  such  materials. 

With  pinpoint  accuracy,  FELs  can  destroy  precise  targets  in  rapid  order,  while  not 
harming  non- military  facilities  and  people.  This  is  a  benefit  when  compared  with  other 
weapons,  which  often  destroy  more  than  the  intended  target.  Thus,  the  damage  to 
surrounding  populations  and  facilities  is  minimized. 

With  an  almost  zero  reaction  time,  FELs  enhance  the  survivability  of  the  user, 
especially  in  cases  when  the  threat  is  unexpected.  This  reaction  time  at  long  range 
disables  the  enemy  threat  before  getting  near.  Current  close-in  weapons,  such  as  the 
Phalanx  gun,  suffer  from  dispersion  as  a  result  of  vibrations.  Moreover,  simulations  and 
demonstrations  have  shown  that  the  probability  of  hitting  the  missile  reduces  as  range 
increases.  The  Phalanx  does  not  destroy  the  missile  at  a  distance  sufficiently  far  enough 
away  to  protect  the  ship  and  the  crew  from  serious  damage  [Ref.  9]. 

The  Rapier,  used  in  the  1982  Falklands  campaign,  the  Hawk,  and  the  Arrow 
systems  are  deployed  worldwide.  They  can  provide  short-range  air  defense  of  airfields 
and  smaller  troop  concentrations.  These  systems  are  more  than  20  years  old,  and  even  the 
famous  Patriot  air  defense  system  is  not  as  good  as  thought  to  be.  “Patriot  wasn't  what  I 
would  call  highly  effective.  We  lost  some  lives  to  ballistic  missiles  during  Desert  Storm. 
Twenty- six  Army  soldiers  were  killed  in  the  barracks  in  Dhahran  as  a  result  as  an  attack 
from  a  Scud  missile”  as  stated  by  Lt.  Gen.  Lester  L.  Lyles  [Ref.  10]. 
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FEL  is  designed  in  such  a  way  that  maintenance  is  minimal.  It  can  be  fired 
instantly  and  no  other  maintenance  is  needed.  Furthermore,  it  could  play  a  role  as  a 
“weapon  of  choice”  in  today’s  U.S.  military  because  it  has  the  ability  to  set  power  levels 
on  a  graduate  scale  and  provide  the  user  with  choices  unavailable  in  kinetic  energy 
weapons  as  well  the  ability  to  minimize  casualties  by  permitting  a  selection  of  non- lethal 
options.  It  has  speed- of- light  engagement  with  a  variety  of  targets  and  has  the  potential  to 
produce  a  range  of  precisely  controlled  effects,  as  well  as  the  potential  for  deep 
magazines  and  low  cost  per  shot.  The  absence  of  costly  rounds  or  defensive  intercept 
missiles  causes  the  cost  per  shot  to  drop  to  ~  $2/missile  (using  lgal  of  fuel),  while  the 
RAM  is  ~$0.7  millions/missile  and  the  PHALANX  is  ~$3k/missile. 

B.  FEL  DESIGN 

1.  Application 

The  application  requires  high  power  in  a  small  footprint.  The  proposed  FEL 
consists  of  six  basic  components:  the  electron  injector,  the  linear  accelerator,  the 
undulator,  the  resonator,  the  refrigerator  and  the  electron  beam  dump.  All  these 
components  can  be  packaged  into  a  box  with  dimensions  of  12m  x  4m  x  4m  as  shown  in 
Figure  3. 


LASER  BEAM 


4  m 


12  m 


Figure  3. 


The  configuration  and  operational  parameters  are  provided  in  Table  1. 


Electron  beam  energy 

4=140  MeV 

Relativistic  Lorentz  factor 

y=  275 

Average  current 

/  =  0.6  Amp 

Accelerator  RF  frequency 

a  =750  MHz 

Accelerator  gradient 

17  MeV/m 

Peak  current 

/  =  800  Amp 

Charge 

q  =  0.8  nC 

Electron  beam  radius 

re  =  0.08  mm 

Electron  beam  density  /pulse 

p=  8xl014  /cm3 

Optical  mode  waist  radius 

w0  =  0.08  mm 

Number  of  undulator  periods 

A  =  20 

Undulator  period 

4=3  cm 

Undulator  length 

L=  60  cm 

Undulator  gap 

g  =  1  cm 

Undulator  parameter 

K  =  2 

Resonator  length 

S  =12  m 

Optical  wavelength 

X,  =1  (im 

Quality  factor 

Q  =  4 

Extraction  efficiency 

T)  «1/4A  =  1.25% 

Energy  spread 

Ay  /y  -2/ N  =0.1% 

Rayleigh  length 

z0  =  1.8cm 

Intensity  on  mirrors 

2 1 0  kW/cnf 

System  dimensions 

12mx4mx4m 

Weight 

60  Tons 

Cost 

$60  million 

Table  1.  MW  FEL  Parameters.  (From  [Ref.  11]) 
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Some  of  the  challenging  issues  are  intensity  in  the  mirrors,  effective  recirculation 
of  the  84  MW  electron  beam  and  beam  transportation.  A  deployable  MW-class  FEL  will 
need  a  sophisticated  feedback,  control  and  diagnostic  system.  Due  to  the  short  turn-on 
time  and  limited  mission  duration,  the  system  must  be  self-diagnosing,  and  self- 
correcting  for  shock  and  vibrationally  induced  misalignments.  [Ref.  12] 

2.  Atmospheric  Propagation 

Although  powerful  at  the  point  of  origin  and  delivered  at  more  than  186,000  miles 
per  second,  a  laser’s  energy  can  be  absorbed  and  diffused  even  in  a  relatively  short 
passage  through  the  atmosphere.  There  are  propagation  issues  that  need  to  be  considered: 

•  Thermal  blooming, 

•  Absorption, 

•  Diffraction, 

•  Scattering. 

It  is  possible  to  overcome  these  difficulties  by  the  use  of  scaling  laws  and  various 
simulations  to  predict  what  will  happen  in  the  MW  class.  These  tools  have  limited 
applications  as  scaling  and  simulations  do  not  always  provide  reliable  results.  The  use  of 
adaptive  optics  is  another  issue  which  may  be  consider  in  our  case. 

a.  Thermal  Blooming 

Random  temperature  variations  caused  by  turbulence  take  place  when  a 
high  power  optical  beam  travels  through  the  atmosphere.  As  a  result,  the  refractive  index 
changes  randomly  and  forms  random  lens  along  the  beam  path.  This  causes  the  optical 
beam  to  spread  as  it  goes  through  the  turbulent  atmosphere.  Thermo  blooming  is  the 
laser- induced  absorption-driven  heating  which  spreads  the  beam  core.  Blooming  is 
nonlinear  and  depends  not  only  on  medium  composition  and  density  but  also  on  beam 
properties  as  well  such  as  wavelength,  pulse  length  and  intensity.  Adverse  weather 
conditions,  humidity,  absence  of  cross  winds  and  the  smoke  from  combat  environments 
can  also  degrade  the  beam  quality. 

b.  Absorption 

One  of  the  major  limitations  of  the  transfer  of  laser  energy  from  the  source 
to  the  target  is  the  absorption  of  the  laser’s  energy  by  the  earth’s  atmosphere.  Light 
passing  through  an  optical  system  can  be  attenuated  by  absorption  as  shown  in  Figure  4. 
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Figure  4.  Light  Absorption. 


The  exponential  law  of  absorption  is  the  basic  working  relationship,  but  specific  terms 
such  as  absorbance  and  transmittance  are  widely  used  as  well.  The  differential  absorption 
can  be  expressed  as 

d<P  =  -a<Pdz  (2.1) 

where  a  is  the  absorption  coefficient  and  which  upon  integration  from  0  to  z  gives  the 
exponential  law  of  absorption: 


0 


e 


-a  z 


(2.2) 


c.  Diffraction 

Diffraction  limits  affects  propagation  over  a  large  distance  and  therefore 

has  to  be  taken  into  account.  For  distances  larger  than  the  Rayleigh  length,  z  >txw2  IX, 
where  w  is  the  radius  of  the  source  aperture,  we  are  able  to  use  the  Fraunhofer 
approximation.  At  a  wavelength  of  1.0  (im  and  an  aperture  width  of  2.5  cm  (1  inch),  this 
distance  is  z  >  1250  meters.  If  we  deal  with  a  circular  aperture,  the  intensity  distribution 
can  be  written  as  [Ref.  13] 
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(2.3) 


3(,WA1V(^/Z)T. 

yhz J  1  kwr / z 

where  A-tiw2  ,  w  is  the  radius  of  the  source  aperture,  J1  is  the  Bessel  function  of  order 
1,  and  r  the  radius  coordinate  in  the  observation  plane.  The  intensity  distribution  in 
Equation  (2.3)  is  referred  to  as  the  Airy  pattern  after  G.  B.  Airy  who  first  derived  it.  Since 
sin0  =  r /  z,  the  irradiance  can  be  written  as  a  function  of  0  ,  and  Equation  (2.3)  gives  the 

intensity  3  at  distance  z  equal  to 

3  =  43(0) 
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Figure  5. 


The  fundamental  mode  (TEMoo),  from  a  circular  aperture,  is  shown  in 
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Figure  5.  Intensity  Pattern  for  Circular  Aperture. 


To  find  the  width  of  the  central  lobe,  measured  along  the  horizontal  axis, 
the  numerator  of  Equation  (2.4)  where  it  is  zero  (1st  of  minima)  must  be  found.  From 
tables,  it  can  be  seen  that  /,  (u)  =  0  when  u  =  3.83  so  that  the  width  of  the  central  lobe  at 
distance  z  is 

D  =  1.22—.  (2.5) 

w 
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Equations  (2.3)  and  (2.5)  strongly  suggest  short  wavelength  for  a  directed  energy  weapon 
as  the  intensity  at  the  target  is  inverse  proportionally  of  the  square  of  the  wavelength, 
while  the  width  where  the  beam  is  focused  is  proportionally  to  the  wavelength. 
Furthermore,  circular  apertures  are  preferable  as  84%  of  light  arrives  within  the  central 
lobe  (Airy  disk),  and  91%  within  the  bounds  of  the  second  dark  ring  [Ref.  6].  For  the 
FEE  weapon  with  X  =  1  pm  wavelengths  laser  light  from  a  w  =  0.2  m  aperture  on  the 
ship,  the  spot  size  at  range  z=  5  km  could  be  as  small  as  wT  =cm.  To  avoid  thermal 
blooming,  to  damage  a  larger  target  area,  the  beam  only  will  be  focused  to  a  10  cm  spot. 
Diffraction  at  X  =  1  pm  is  not  a  problem. 
d.  Scattering 

Rayleigh  scattering  refers  to  the  scattering  of  hght  off  molecules  in  the  air, 
and  from  particles  up  to  about  a  tenth  of  the  wavelength  of  light.  The  strong  wavelength 
dependence  of  Rayleigh  scattering  favors  the  short  wavelength  X  since  the  scattered 

intensity  I  is  proportional  to  X  1  [Ref.  14].  For  example  the  Rayleigh  scattering  at  400  nm 
is  40  times  greater  than  at  1pm  for  equal  incident  intensity. 

For  particles  sizes  larger  than  a  wavelength,  Mie  scattering  predominates. 
This  scattering  produces  a  pattern  like  an  antenna  lobe,  with  a  sharper,  more  intense 
forward  lobe  for  larger  particles.  Mie  scattering  is  not  strongly  wavelength  dependent. 
Both  scattering  types  are  presented  schematically  in  Figure  6. 


COMBAT  ENVIRONMENT 


Rayleigh  Scattering  MieScattering 


Mie  Scattering, 
larger  particles 

> 

-► 

*- 

> 


Incident  light 


Figure  6.  Scattering  from  Air  Molecules. 
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The  atmospheric  transmittance  over  a  wavelength  range  is  a  complex 
problem  and  extremely  crucial  to  directed  energy  because  it  determines  the  choice  of  a 
suitable  wavelength.  Absorption- driven  thermal  blooming,  scattering  and  turbulence  can 
be  major  limitations  to  performance  and  can  lead  to  the  dispersion  of  the  beam  power. 

The  atmospheric  transmittance  problem  has  been  described  by  computer 
modeling  codes  such  as  FASCODE,  LOWTRAN,  MODRAN  and  MOLLY.  All  the 
above  codes  allow  the  user  to  insert  weather  conditions  and  to  select  from  several 
different  environments.  MOLLY  is  a  time- dependent  computer  simulation  of  adaptively 
compensated  laser  propagation  through  turbulence  and  thermal  blooming.  A  graph  of 
typical  atmospheric  transmittance  versus  wavelength  using  this  simulation  is  shown  in 
Ligure  7.  Notice  there  are  specific  windows  where  the  attenuation  is  much  less,  and 
therefore  it  is  preferable  to  operate  the  laser  in  those  regions. 


Maritime  Atmosphere 
Marine  Aerosol,  Horizontal  Path 


Ligure  7.  Atmospheric  Transmittance  (Lrom  [Ref.  15]). 

As  can  be  seen  from  the  above  figure,  both  COIL  and  MIRACL  lasers 
perform  poorly  in  a  maritime  environment  and  suffer  from  absorption.  Wavelengths 
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X  =  1.06  [tm,  X  =1.25  (im,  X  =1.62  (im  are  attractive  for  the  MW-class  because  as  the 
absorption  coefficient  is  small,  a  <  0.001km"1.  For  the  value  of  a  =  0.001km"1  we  have 
almost  no  reduction  of  the  initial  power,  see  Equation  (2.2),  and  d>0  ~<f>F .  Note  that  the 
width  of  each  window  is  relatively  small  (~0.2|im),  and  a  small  deviation  causes 
excessive  dissipation  and  spreading  due  to  atmospheric  conditions,  which  results  in 
degraded  range  and  beam  quality. 

3.  Adaptive  Optics 

Current  acquisition  and  pointing  systems  are  capable  of  tracking  maneuverable 
supersonic  targets  and  holding  the  optical  beam  to  a  small  area  of- 100cm  .  Adaptive 
optics  can  be  used  to  transfer  the  lethal  beam  and  focus  it  in.  In  addition,  adaptive  optics 
are  used  to  overcome  propagation  affected  by  turbulence..  The  advantage  of  using 
adaptive  optics  is  shown  in  the  simulations  results  in  Figure  8. 


no  adaptive  optics  adaptive  optics 

Figure  8.  Function  of  Adaptive  Optics  (from  [Ref.  15]). 

On  the  left  hand  side,  no  adaptive  optics  are  used  and  thermal  distortions  in  the 
atmosphere  spread  the  laser  beam.  On  the  right  hand  side  the  improvement  is  obvious 
after  the  application  of  adaptive  optic  mirror  at  source.  Thermal  distortions  and  moderate 
thermo  blooming  are  largely  corrected  and  allow  the  delivery  of  the  lethal  beam  to  the 
target.  The  system  performance  is  degraded  if  the  thermal  distortions  and  the  thermo 
blooming  effect  are  more  than  moderate.  In  the  FEE  application  discussed  here, 
turbulence  is  much  less  than  shown  in  Figure  7  so  that  adaptive  optics  may  not  be 
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required.  If  it  is  required  the  adaptive  optics  adds  a  cost  of  $3M  to  the  entire  FEL  system 
[Ref.  8], 

C.  FEL  SAFETY 

1.  Radiation  Concerns 

Radiation  shielding  is  required  for  the  electron  beam  dump.  The  FEL  with  the 
energy  recovery  system  allows  the  electron  beam  to  reenter  the  RF  cavities  180  degrees 
out  of  phase  with  the  accelerating  fields.  This  energy  is  recovered  as  RF  power  and  used 
to  accelerate  subsequent  electron  pulses  to  repeat  the  cycle.  The  residual  energy  after 
deceleration  is  about  7  MeV  [Ref.  16]  and  is  below  the  threshold  of  generating  neutron 
radiation  when  it  is  dissipated  in  the  beam  dump.  Only  beryllium  (Be)  and  Osmium  (Os) 
have  a  neutron  production  threshold  below  7  MeV  and  therefore,  as  long  as  beryllium 
and  osmium  are  excluded  from  beam  dump  structural  materials,  neutron  generation  is  not 
an  issue  [Ref.  17].  Without  the  production  of  neutrons,  much  less  shielding  is  required. 
Other  human  safety  issues  include  the  need  for  eye  protection  when  personnel  may  be 
topside  or  have  visual  access  to  the  topside  scattering  off  the  sea  surface  during  an 
engagement.  In  such  cases,  safety  goggles  may  be  required  [Ref.  18]. 

2.  Helium  Concerns 

The  cryogenic  systems  require  the  use  of  liquid  helium.  Helium  is  non¬ 
flammable,  non-corrosive,  non-toxic  and  not  listed  as  a  marine  pollutant  by  Department 
of  Transportation  (DOT).  When  the  FEL  is  running  properly,  the  helium  is  in  a  closed 
loop  refrigerator  and  totally  contained.  Thus,  there  are  no  safety  concerns.  However, 
because  the  helium  is  cryogenic,  a  vacuum  is  used  to  insulate  the  piping  and  storage 
containers  (dewars).  If  a  major  vacuum  leak  should  occur,  the  helium  boils  off  quite 
rapidly  leading  to  safety  concerns: 

•  It  can  over  pressure  the  storage  tank  or  piping  since  helium  expands  750 
times  when  going  from  a  room  temperature  liquid  to  a  gas.  If  a  suitable 
way  is  not  provided  for  the  gas  to  escape,  the  vessel  could  explode. 

•  If  the  expanding  gas  should  come  into  contact  with  someone,  cryo-bums 
can  occur  because  it  is  very  cold  when  first  released. 

•  If  the  helium  goes  into  a  confined  space,  it  can  displace  the  oxygen  and 
lead  to  suffocation  of  personnel  in  those  compartments.  If  air  with  less 
than  12%  oxygen  is  breathed,  unconsciousness  can  occur  in  9  seconds 
with  brain  damage  occurring  in  5  minutes.  If  the  level  drops  to  8%,  death 
is  imminent  in  minutes. 
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The  precautions  are: 

•  Provide  burst  disks  on  the  system  connected  to  a  vent  so  the  helium  can 
safely  escape. 

•  In  the  case  of  an  expanding  plume,  do  not  stick  a  body  part  into  it. 
Protective  clothing  is  required. 

•  Provide  a  means  for  the  helium  to  escape  if  it  wants  to  rise.  If  there  is  a 
pipe  in  the  ceiling,  it  will  escape  up  the  pipe  and  disperse.  To  prevent 
leakage  into  other  compartments  have  a  ledge  at  the  top  of  the  door  below 
the  ceiling  level  to  prevent  the  flow.  Also,  oxygen  monitors  on  the  ceiling 
to  send  out  an  alarm  if  oxygen  levels  drop  below  18%  are  required. 

The  total  volume  of  liquid  helium  in  the  linear  accelerator  will  be  less  than  5000 
liters,  making  approximately  5  megaliters  of  gas  available  which  could  easily  fill  a  room 
with  155m  volume.  Nevertheless,  helium  is  much  safer  than  liquid  nitrogen.  It  goes  to 
the  ceiling  and  escapes.  Nitrogen  goes  to  the  floor  and  cannot  be  seen  once  it  warms  [Ref. 
19].  Helium  vents  would  be  included  in  the  weapon  design. 

D.  CONCLUSIONS 

1.  Challenges 

Theoretical  calculations  have  shown  that  FELs  have  the  potential  to  scale  to 
megawatts,  but  it  is  an  engineering  challenge  to  improve  reliability  and  reduce  the  size 
and  the  cost  of  an  FEL.  For  instance,  a  critical  concern  is  improving  electron-beam 
brightness  as  beam  current  is  increased.  Additionally,  new  functions  must  be  introduced 
into  FEL  technology  to  advance  the  MW  system.  The  radius  of  curvature  of  the 
resonators  mirrors  and  a  cooled  sapphire  window  must  be  controlled  to  tight  tolerances. 
Furthermore,  atmospheric  propagation  effects,  and  the  thresholds  at  which  these  effects 
are  significant,  must  be  well  understood.  This  understanding  will,  to  a  large  extent,  define 
the  required  laser  weapon  system.  Performance,  downsizing  and  cost  are  also  crucial 
considerations  and  are  illustrated  in  Figure  9. 
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HEL  Systems  Considerations 


Figure  9.  Acceptable  System  Diagram  (After  [Ref.  20]). 

After  the  workshop  in  Virginia  in  June  2001,  the  synopsis  is  that  a  MW-class 
infrared  free-electron  laser  (IRFEL)  appears  to  be  a  challenging,  but  feasible,  proposition 
[Ref.  21]. 

2.  Resources  Required 

The  DoD  investment  strategy  is  to  fund  those  technologies  that  have  the  potential 
to  penetrate  and  bypass  identified  technology  barriers.  In  practice,  in  the  face  of  funding 
pressures,  inadequate  funding  is  being  provided  to  a  wide  variety  of  programs.  These 
programs  should  be  replaced  with  focused,  sequential  developments  funded  at  the  level 
of  effort  needed  to  create  real  progress.  It  is  also  believed  that  considerably  more  funding 
on  the  level  of  $100-150  million  per  year  is  needed  [Ref.  22]. 

Currently,  the  money  spent  in  this  area  in  FY01  was  $30  million,  while  the 
funding  for  the  DoD  high-energy  laser  programs  was  $474  million  [Ref.  20].  Studies  at 
TJNAF  and  simulations  presented  in  Chapter  IV  have  indicated  that  scaling  the  system  to 
100  kW  of  power  is  feasible  and  straightforward.  Funding  is  already  being  supplied  by 
the  Office  of  Naval  Research  (ONR)  and  additional  funding  is  expected  from  the  U.S.  Air 
Force. 
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in.  THE  PHYSICS  BEHIND  FEL 


A.  INTRODUCTION 

A  free  electron  traveling  undisturbed  will  not  spontaneously  emit  energy. 
However,  if  the  electron  accelerates,  it  must  absorb  energy.  If  it  decelerates,  it  must 
radiate  energy.  For  example,  if  an  electron  increases  its  velocity  in  a  particle  accelerator, 
it  does  so  by  absorbing  energy  from  an  RF  electromagnetic  field  supplied  by  the 
accelerator.  Its  kinetic  energy  is 

Ke=(y  - X)mec2  (3.1) 

where  y  is  the  Lorentz  factor  and  me  is  the  electron  rest  mass.  The  term  mec2  is  called  rest 
energy  (0.511  MeV).  The  Lorentz  factory  is  given  by 

y=-7  l—  0.2) 

Vi -PP 

where  p  =v  /  c  is  the  dimensionless  electron  velocity. 

In  an  FEL,  the  electron  trajectory  is  directed  with  bending  magnets  into  the 
undulator.  Along  the  path  of  the  undulator,  the  electron  starts  to  wiggle  in  the  transverse 
direction  as  soon  as  it  passes  through  each  alternating  magnetic  field.  During  the 
wiggling,  the  electron  emits  radiation  into  a  narrow  cone  of  opening  angle  "0  [Ref.  23], 
where 

=y_1  (3.3) 

This  type  of  radiation  should  be  distinguished  from  Bremsstrahlung  or  broad  band 
synchrotron  radiation  which  usually  consists  of  several  harmonics  in  addition  to  the 
fundamental  [Ref.  24].  This  thesis  deals  primarily  with  a  pulsed  FEL  operated  as  an 
oscillator.  A  radio- frequency  (RF)  electron- gun  is  used  to  inject  electron  pulses  into  the 
accelerator.  The  radiation  is  stored  in  a  cavity  and  amplified  over  many  passes.  Other 
types  of  FELs  either  use  continuous  wave  (CW)  power  or  amplify  the  radiation  in  a 
single  pass  and  are  called  amplifiers  [Ref.  25].  Table  1  lists  not  only  the  different  types, 

but  also  the  design  parameters  for  existing  and  proposed  FELs  [Ref.  26]. 
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Relativistic  Short  Wavelength  Tree  Electron  Lasers  (2000) 
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RF  -  RF  Linac  Accelerator 
SR  -  Electron  Storage  Ring 
A  -  III  Amplifier  " 

SP  -  Smith-Purcell  Oscillator 


MA  -  Microtron  Accelerator 
EA  -  Electrostatic  Accelerator 
O  -  EEL  Oscillator 
S  -  SASE  FEL 


Table  2.  The  Short  Wavelength  FEL  in  2000  (From  [Ref.  26]). 
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The  first  column  lists  the  operating  optical  wavelength  X  which  is  given  by  [Ref.  25] 


X  = 


K0- +K2) 

2y2 


(3.4) 


where  X()  is  the  undulator  wavelength,  and  K  is  the  undulator  parameter  listed  in  the 
following  column.  The  relativistic  Doppler  shift  and  Lorentz  contraction  due  to  the 
electron  motion  is  included  in  the  calculation  of  the  radiated  optical  wavelength  in 
Equation  (3.4).  This  is  called  the  resonance  condition  and  occurs  when  the  electron  slips 
behind  the  optical  wave  by  one  optical  wavelength  as  the  electron  traverses  an  undulator 
wavelength.  From  the  simple  relationship  of  (3.4),  it  is  seen  that  an  FEL  can  easily  be 
designed  to  operate  at  a  broad  span  of  wavelengths,  including  wavelength  regions  where 
there  are  no  powerful  light  sources.  For  example,  no  conventional  laser  operates  in 
ultraviolet  or  an  extreme  ultraviolet  range  of  the  electromagnetic  spectrum  (300  to  30 
nm)  and  no  conventional  laser  operate  in  the  far  infrared  (1mm  to  20  (im)  [Ref.  27]. 


The  second  column  lists  the  normalized  electron  pulse  length 


(3.5) 


where  i  e  is  the  electron  pulse  length  and  c  the  speed  of  light.  From  this  column,  it  can 
clearly  be  seen  that  the  majority  of  the  FEFs  create  micropulses  and  only  a  few  work  CW 

The  third  column  provides  the  electron  beam  energy  E  as  it  comes  from  the 
accelerator,  ranging  from  a  few  Me  Vs  up  to  1  GeV  in  ELETTRA  while  the  proposed 
electron  beam  energy  is  up  to  30  GeV  in  the  TESLA  FEE  facility.  The  next  column  shows 
the  average  peak  current  /  ranging  from  a  few  milliamperes  to  300A  for  the  existing 
FELs  and  up  to  5000A  for  the  proposed  FELs. 

The  fifth  and  sixth  columns  are  parameters  related  to  the  undulator  design.  The 
number  of  undulator  periods  is  N,  and  the  undulator  wavelength  is 


(3.6) 
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where  L  is  the  undulator  length.  Note  that  in  a  few  cases,  in  the  iV-column,  there  are 
multiple  undulator  sections,  which  are  referred  to  in  the  FEL  klystron  undulator  designs 
[Ref.  26] .  The  next  column  also  depends  on  the  design  and  is  the  undulator  parameter 


K  = 


eBX0 
2n  mec 2 


(3.7) 


where  e  is  the  electron  charge  magnitude  and  B  is  the  rms  undulator  field  strength.  Take 
special  note  that  the  values  of  K  ~  1 .  For  K  » 1 ,  there  will  be  radiation  in  many  higher 
harmonics  [Ref.  24],  and  for  K«l,  FEL  gain  would  be  small. 

The  last  column  lists  the  accelerator  type.  Note  that  the  majority  are  RF  linear 
accelerators  and  oscillators,  as  in  the  experiment  and  simulations  described  in  the 
chapters  that  follow. 

B.  FEL  INTERACTION-PENDULUM  EQUATION 

In  this  thesis,  as  well  as  in  most  cases,  Coulomb  forces  between  the  electrons  are 
small  because  the  beam  energy  is  large.  Relativistic  electrons  streaming  in  the  z  direction 
and  interact  with  a  helical  undulator  magnetic  field, 

B  =  B 0  (cosk^sinA^O)  (3.8) 


where  B  is  the  magnetic  field  of  the  undulator  and  k0  =  27t  /  \  is  the  undulator 
wavenumber.  The  electron  Lorentz  force  equation  is  [28] 

£®)=_JL(£  +  |Jxs)  (3.9) 

dt  mec ' 

Initially,  we  find  the  electron  trajectory  through  the  undulator  by  assuming  that  there  is 
no  radiation  E  =  0.  The  electron  energy  change  is  given  by 

=  (3.10) 

dt  mec 


which  means  that  d y  /  dt  is  zero  according  Equation  (3.10)  for  the  E  =  0  case  ,  so  that  y 
would  be  constant.  The  term  (3  x  B  in  Equation  (3.9)  is 

P  x  B  =  sin  k0z )  +  y(B0  p,  cosk0z)  +  z(B0 Px  sink0z~  50Py  cosk0z)  (3.1 1) 
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Then,  Equation  (3.9)  can  be  separated  into  components  and  becomes 


rfP,  _  eB$z  • 


dt  ymec 


sin  k()z, 


(3.12) 


which  by  integration  gives 


P  v  = e-\r cos  Kz 

Y mec  k0 


(3.13) 


since  dz  / dt  =  P_c,  and  where  it  is  assumed  that  the  ejection  of  electrons  is  perfect 
and  the  constant  of  integration  is  zero.  Substitute  k0  =2n  /\  Equation  (3.13)  becomes 


Pa  =  “  z~~^~Tcos  k0Z 

27ty  mec 


(3.14) 


Define  the  undulator  parameter  as  K  =  eB0X0  /  2tv/  me2  so  that  Equation  (3.14)  becomes 


PA  =  -— COS^qZ. 

y 


(3.15) 


Similarly  for  the  y- component,  we  have 


Pv  =-— (sin  k0z). 


(3.16) 


The  vector  addition  of  the  components  in  Equations  (3.15)  and  (3.16)  gives  the  electron 

— ► 

velocity  in  the  transverse  direction  [1 L  =  ( (3^ ,  (3  ,0)  .  Thus, 


P±  = - (cosfc0z,sin  kQz,0) 

y 


(3.17) 


Now  that  we  have  the  electron  trajectory,  we  add  radiation  to  find  the  microscopic 
motion.  Taking  into  account  the  radiation  created  in  the  undulator,  consider  a  circularly 
polarized  plane  wave  for  the  optical  field  of  the  form 


Es  =  E(cos\\r ,  -sin\|/,0), 


(3.18) 
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Bs  =  ii(sin\|/,cos\|/,0). 


(3.19) 


Here  \\i  =  kz  -tot  -ftp  ,  E  is  the  electric  and  magnetic  field  amplitude  in  cgs  units,  k  =  to/c 
is  the  optical  wavenumber,  and  cp  is  the  optical  phase. 

With  radiation  included  Equation  (3.9)  becomes 

—  = - — £(Pxcos\|/-Pvsin\|/).  (3.20) 

dt  mec 

Substitution  of  Equations  (3.15)  and  (3.16)  into  (3.20)  yields 

pKF 

y= - cos£  +(p),  (3.21) 

y  mec 


where  y  =dy  / dt  and  'Q  +cp  =\|/  +k0z  =  (k+k0)z-(Ot  +  q> .  C,  is  defined  as  an  electron 
phase,  which  describes  the  electron  position  in  the  combined  undulator  and  optical  fields 
( /c0and  k  terms).  Note  that  in  Equation  (3.21),  the  rate  of  change  in  energy,  y ,  is 
proportional  to  the  electric  field  E  and  the  undulator  parameter  K,  but  inversely 
proportional  to  the  energy  of  the  electrons  y.  As  the  energy  of  the  incoming  electrons  is 
increased,  the  coupling  with  the  optical  field  is  reduced. 


Expanding  Equation  (3.2),  the  Lorentz  factor  reads 

y-2=l-P;-(3;, 

and  substituting  Equation  (3.17)  in  Equation  (3.22)  yields 

y-2(l+Z2)  =  l-p2. 

Taking  the  time  derivative  on  both  sides  of  Equation  (3.23)  gives 

-2Y-3Y(l+jrJ)  =  -2PJ!=A  =  rflA 

y  1 +K 

Recall  that  the  electron  phase  is 

C  =  (k  +kQ)z  -to t. 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


and  by  differentiating  Equation  (3.25)  twice  gives 


26 


C  =  [(k + ko)z]  -  =  [(k  +  K )  Pz  ] _(0  ’ 

and  by  inserting  the  result  of  (3,  in  Equation  (3.24)  we  get 

y  _  yjX 

y  (1  +  K~)(k  +k0)c 


(3.26) 

(3.27) 


(3.28) 


For  relativistic  electron,  k  »  k0  and  (1  ~  1 ,  so  that  Equation  (3.28)  can  be  written  as 

(3-29) 

y  (1 +K~)kc  {l+K-)0) 

The  resonance  condition  of  Equation  (3.4),  with  X  is  replaced  by  271  c/ co  can  be  written 
as 


co  = 


co()2y 

1 +  K2  ' 


(3.30) 


Substituting  Equations  (3.21)  and  (3.30)  in  Equation  (3.29)  and  solving  for  £  gives 

••  2m  eKF 

C  =  —^—§-cos((:  +tp).  (3.31) 

mecy 

Dimensionless  parameters  are  introduced  to  make  the  model  simpler  and  more 
physically  meaningful.  The  dimensionless  time  x  =ct  /  L  describes  the  electron  and 
optical  evolutions  that  take  place  as  X  ranges  from  0  — >  1 .  The  time  derivatives  now  are 

(..)  =  d(..)/dx  instead  of  (..)  =  d(..)/ dt .  The  dimensionless  optical  complex  field  is 
defined  as 

a=\a\ei(f  (3.32) 

II  2  2 

where  a  =4nNeKLE /y  me  is  the  dimensionless  optical  field  amplitude.  Applying  all 

the  new  dimensionless  quantities,  Equation  (3.31)  takes  the  form  of  a  simple  pendulum 
equation 
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(3.33) 


£  =  |a|cos(^  +cp). 

Equation  (3.33)  governs  the  electron’s  microscopic  phase  dynamics  relative  to  the 
optical  field  inside  the  undulator.  Therefore,  the  complicated  motion  of  the  electrons  is 
described  by  this  simple  and  well-known  pendulum  equation  The  positive  or  negative 
acceleration  in  Equation  (3.33),  is  proportional  to  values  of  the  cosine  term  and  results  in 
the  bunching  of  the  electrons  within  an  optical  wavelength.  Consider  a  group  of  electrons 
uniformly  distributed  in  phase.  The  electrons  with  phases  in  the  range 
-Til 2  <C,  +cp  <n  / 2  will  gain  energy  and  speed  up,  while  the  electrons  in  the  range 

7t/2<^+(p<37t/2 will  experience  a  loss  in  energy  and  slow  down.  As  a  result,  the 
electrons  have  the  tendency  to  bunch  together  near  C,  =tz  /  2  .  This  bunching  coupled  with 
the  wave  equation  contributes  to  a  net  energy  transfer  to  the  optical  field. 

The  dimensionless  field  amplitude  |a|  of  Equation  (3.33)  determines  the  electron 
bunching  rate  along  the  undulator.  If  |a| » 7t ,  the  optical  field  is  strong  and  the  bunching 
occurs  quickly.  If  |a|«7t  ,  the  optical  field  is  weak  and  the  change  in  the  electron  phase 
velocity  is  small  [Ref.  24]. 

The  wave  equation  that  drives  the  complex  optical  dimensionless  field  of 
Equation  (3.32)  can  be  put  into  the  following  simple  form  [Ref.  25] 

a  =  -j  <  e-K  >  (3.34) 

where  the  dimensionless  electron  beam  current  j  =  8N(eizKL )2  p  /y  3mec2 ,  p  is  the  actual 
electron  beam  density,  and  ^..)  in  the  wave  equation  indicates  an  average  over  all  the 
electrons.  This  average  will  be  non- zero  only  when  the  electrons  are  bunched.  The  rate  of 
change  of  the  optical  field  is  proportional  to  the  dimensionless  current  j  . 

The  pendulum  Equation  (3.33)  and  the  wave  Equation  (3.34)  are  coupled  together 
through  the  dimensionless  current  j.  If  j  <  1 ,  the  gain  is  low,  and  if  j  » 1  the  gain  is 
high.  Both  the  pendulum  and  wave  equations,  even  in  this  simple  form,  are  valid  in  high 
and  low  gain  regimes,  as  well  as  in  weak  and  strong  optical  fields. 
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C.  PHASE  SPACE  ANALYSIS 

The  electron  evolution  and  the  optical  wave  are  coupled  accordingly  in  Equations 
(3.33)  and  (3.34).  The  phase- space  evolution  of  sampled  electrons  is  used  to  describe  the 
motion  of  the  electrons  over  an  optical  wavelength.  Each  electron  is  started  with  its  initial 
conditions  =  £  (0)  and  v0  =V  (0)  at  the  beginning  of  the  undulator  (x  =  0) .  The  phase 
velocity  is  given  in  dimensionless  notation  as 

=L[(k+k0)$z-k].  (3.35) 

dx 

At  resonance,  V=0,  so  solving  (3.35)  for  (1  and  using  the  resonance  condition 
\  =\(\  +  K2)  1 2y2  and  the  approximation  that  y  » 1 : 


k 

k  +At() 


1 +K2 


(3.36) 


Substituting  (Y  into  Equation  (3.35)  yields 


V  =  L 


1 +K2 

(k+kQ)(l-—^)-k 

2y 


f\  +  K2  A 

j 


(3.37) 


using  k0«  k.  Now  differentiating  Equation  (3.37)  with  respect  to  y,  substituting 
L  =  N\  ,  and  applying  the  resonance  condition  gain  gives: 


Av  =4t zN 


(3.38) 


U  J 

From  Equation  (3.38),  it  is  clear  that  the  change  in  phase  velocity  is  related  to  the 
change  in  the  electron  energy.  Figure  10  shows  the  phase  space  evolution  for  electrons  at 
resonance  and  above  resonance.  In  the  first  case,  electrons  are  injected  into  the  undulator 
with  initial  electron  velocity  v  0  =  0  .  They  appear  as  yellow  spots  (evolution  denoted  by 
yellow  going  to  red).  Note  that  the  electrons  between  -n/2  <£  < 7t  / 2  increase  their 
phase  velocity  by  absorbing  energy  from  the  optical  field.  An  equal  number  of  electrons 
between  n  /  2  <C,  <3%  /  2  decelerate  and  transfer  energy  to  the  optical  beam.  Thus,  the 
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net  energy  transfer  is  zero.  On  the  top  right  of  the  figure,  where  the  output  of  the  optical 
gain  G(t  )  is  located,  it  can  be  clearly  seen  that  a  resonant  electron  beam  has  no  gain. 


****  p£|_  p|-,ase  space  Evolution  **** 
j=l  a  o=7r  vo=0  N  =  36 


(a) 


****  PEL  Phase  Space  Evolution  **** 
j=l  a 0=7t  v0=2.6  N  =  36 


1 2 


GJa.in 


□  .3 


(b) 


Figure  10. 


The  Phase-Space  Evolution  at  Resonance  (a),  at  Optimum  Initial  Velocity  (b). 
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On  the  other  hand,  the  bunching  that  occurs  near  the  relative  phase  C,  +cp  =rc  /2, 
drives  the  optical  phase  cp  in  Equation  (3.33).  The  optical  phase  cp(x)  is  plotted  below 
the  gain  and  it  is  relatively  large,  (-0.1). 

Electrons  injected  slightly  above  resonance  at  v0  =  2.6,  are  shown  in  the  bottom 

part  of  Figure  10.  In  the  phase  space  evolution  for  x  =  0  — » 1  along  the  undulator, 
bunching  near  C,  =7C  occurs.  This  bunching  indicates  that  a  significantly  larger  number 
of  electrons  transfer  energy  to  the  optical  field  and  the  gain  is  increased  to  G  =12.5%  . 

The  curve  is  shown  in  red  in  the  same  figure  is  called  the  “separatrix”.  It  separates 
the  closed  and  open  orbits  in  phase  space  (phase- space  paths)  and  connects  unstable  fixed 
points  at  C,  = -K  12,  371/2  and  v  =0  .  The  separatrix  points  ( C  v ,V  s )  are  given  by  [Ref. 
25] 

v;  =  2|«|[1  +sin(C,  +cp)].  (3.39) 

I  |l/2 

The  peak-to-peak  height  of  the  separatrix  is  4  a  .  Electrons  outside  the 

separatrix  follow  open  orbits;  those  inside  the  separatrix  are  trapped  in  closed  orbits. 

D.  GAIN 

From  the  wave  Equation  (3.34),  it  is  clear  that  if  j  =  0  ,  then  a  -  0  and  there  is  no 
gain.  If  the  dimensionless  current  is  small  j  «  n  ,  then  a  small  change  in  the  optical  field 
amplitude  occurs,  and  the  gain  is  low  since  the  electron  phase  velocity  does  not  change 
much  during  the  evolution  through  the  undulator.  Today,  the  majority  of  the  operating 
FELs  are  low-gain,  low-current  FELs,  and  based  on  the  following  assumptions,  it  is 
possible  to  derive  the  weak  field  gain. 

•  Assume  weak  fields:  a«  n,  cp=0 

•  Assume  low  gain:  a  ~  0=^>  a  =  a0,  cp  ~  0 

With  these  assumptions,  the  pendulum  equation  can  be  solved  analytically  by  expanding 
C,  and  v  in  a  power  series. 

The  zero  order  expansion  is 
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(3.40) 


C(0)=Co+V0T 

V(0)  =V0. 


To  first  order  in  a0 ,  where  a0  =  a(0)  ,  the  pendulum  equation  gives 


C  U)  =-^t[cos(Co  +V0x  )-cosC0  +V0x  sin^0] 
v  01  =[sin(C0  +v0x  )-sin£0]. 


(3-41) 


In  order  to  have  net  energy  transfer  from  the  electrons  to  the  optical  beam,  we 
must  satisfy  (v)-v0  ^0 ,  where  (v)  is  the  average  phase  velocity  over  all  the  electrons. 

The  average  of  the  (v  (1)^  in  Equation  (3.41)  is  zero  because  the  average  of  the  sin 


function  is  (sin)  =  0  which  means  that  just  as  many  electrons  gain  energy  as  lose  energy. 
Therefore,  second  orders  terms  are  needed  to  obtain  non-  zero  gain,  so  that  we  find 


V  (2 >  =  ■%  (cos(2^  0  +  2v  0x )  -  cos(2£0 ) )  +  cos(v  0x )  - 1  -  v  0x  sin(£  0)cos(£  0  +v  0x ) 

v0  L  4 


1 


(3.42) 


The  average  of  the  second  order  term  is 


=v0+- 


V 


0  L 


cos(v  0x )  - 1  +  0x  si  n  (v  ()x ) 


(3.43) 


The  above  dynamics  of  the  electron  position  is  related  to  the  energy  transfer,  with  the 
relation  Av  =4n A/' (Ay  /y) .  Both  sides  are  averaged  to  find  the  average  change  in  the 
electron  beam  energy 


Ay  me1  =ymc 


((v)-v.) 

47 iN 


(3.44) 


Energy  conservation  is  used  to  calculate  the  gain  G(x ) .  Gain  is  defined  as  the  ratio  of  the 
energy  transfer  from  electrons,  Equation  (3.44),  divided  by  radiation  beam  energy  in  cgs 
units,  2E]dV  1  8ji  ,  in  a  volume  dV,  where  Es  is  the  electric  optical  field.  The  filling 
factor  F  is  defined  as  the  ratio  of  the  electron  beam  area  divided  by  the  optical  mode  area 
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(3.45) 


where  rb  is  the  electron  beam  radius  and  w0  is  the  optical  beam  radius  at  the  waist.  The 
number  of  electrons  in  a  volume  dV  is  given  by  p  FdV ,  where  p  is  the  electron  density. 
Thus  the  gain  at  time  x  is  given  by 

(p FdV )y me2  ((v)-Vn)/47CiV 

G(t  )  =  -  — - - — - -1 - .  (3.46) 

2E-dV/Sn 

In  the  exchange  of  the  energy,  the  second  order  contribution  of  Equation  (3.43)  is  used. 
Also,  the  definition  of  dimensionless  current  j  =  8N(&lKL)2 p  F /y3mc2  is  used.  Using 
these  relations  in  to  Equation  (3.46)  yields 


G(x)  =^-{2-2cos(v0x)-v0x  sin(v0x)}.  (3.47) 

v0 

From  the  above  equation,  the  strong  relationship  between  the  gain  of  the  optical  field  and 
the  initial  phase  velocity  v0  can  be  seen. 


Figure  11  shows  the  output  of  simulations  in  the  weak- field,  weak-current  regime. 
The  final  gain  G(x  =  1) ,  as  well  as  the  final  optical  phase  shift  cp(x  =1)  versus  the  initial 

phase  velocity,  v0,  is  plotted. 


33 


0.1324 


****  Gain  and  Phase  Curves  **** 

j=l  a  =1  o=0 

J  o 

N=4 1  8Q-0 


0.0 

-0.1324 

0.079 

0.0 


-0.079 


Figure  11.  Weak-Field,  Low-Current  Gain  G(V 0)  and  Optical  Phase  Shift  cp(v0)  Spectra. 

The  simulation  result  is  in  total  agreement  with  the  analytic  results  of  Equation 
(3.47).  The  initial  phase  velocity,  which  is  a  function  of  the  electron  beam  energy  (recall 
that  Ay  /y  =  Av0  /4n N ),  determines  the  final  gain.  At  the  resonance  condition,  (v0  =  0) , 

the  gain  is  reduced  to  zero.  On  the  other  hand,  the  optical  field  drives  the  optical  phase 
shift,  with  a  peak  value  of  Acp  =0.8  j .  There  is  a  peak  gain  of  G  =  0.13  j  when  the  initial 

velocity  is  v0  =  2.6  and  the  optical  phase  shift  is  reduced  to  Acp  =0.02.  As  the  gain 
spectrum  is  anti- symmetric  around  v0  =  0,  there  are  losses  for  negative  initial  phase 
velocities  (v0  <0)  and  accordingly,  the  maximum  absorption  occurs  at  v0  =-2.6  . 

E.  SHORT  PULSES 

From  the  last  column  of  Table  1  it  is  clear  that  the  majority  of  FELs  use  radio- 
frequency  accelerators  (RF).  Electrons  are  ejected  in  short  pulses  from  an  electron-gun 
with  a  frequency  /  and  are  accelerated.  The  separation  distance  between  the  pulses  is  c/f. 
Thus,  for  /=  750  MHz,  the  pulses  are  separated  by  40  cm.  The  pulse  length  that  is 
currently  used  for  FEL  oscillators  is  often  picoseconds  and  is  comparable  to  the  slippage 

distance  NA,  where  N  is  the  number  of  undulator  periods. 
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In  the  picosecond  scale,  the  short-pulses  effect  the  FEL  interaction.  As  electron 
“buckets”  enter  an  undulator  length  L ,  short  optical  pulses  are  produced  due  to 
spontaneous  emission.  These  optical  pulses  bounce  between  the  mirrors  of  the  resonator 
separated  by  distance  S.  The  optical  pulse  after  a  time  2 Sic  arrives  back  at  the 
beginning  of  the  undulator.  If  the  incoming  electron  buckets  enter  the  undulator  at  the 
same  time  and  coincide  with  the  optical  pulse,  exact  synchronism  occurs.  The 
displacement  between  optical  and  electron  pulses  is  called  desynchronism,  d,  and  in  exact 
synchronism  d  =  0 .  As  the  optical  pulses  enter  the  undulator  simultaneously  with 
electron  pulses, gain  development  is  delayed.  This  effect  is  called  lethargy  [Ref.  25].  The 
electron  pulses  travel  slightly  slower  than  the  optical  pulses.  The  electron  pulses  fall 
behind  and  preferentially  amplify  the  trailing  part  of  the  optical  pulse;  so  the  centroid  of 
the  optical  pulse  travels  slower  than  the  speed  of  light.  After  each  pass,  the  optical  pulse 
falls  farther  behind  the  electron  pulse  and  consequently,  after  some  passes,  the  power 
starts  to  decay. 

Figure  12  shows  the  short  pulse  evolution  at  exact  synchronism,  d  =0  .  The  lower 
left  window  shows  the  electron  pulse  position  for  x  =  0  and  x  =1.  The  electron  pulse 
shape  is  parabolic  and  has  the  form  j(z)  =  j(l-2z2/o2)  for  j(z)  >  0  and  is  zero 

otherwise.  The  normalized  pulse  length  is  o ,  and  the  dimensionless  current  j  =  5  .  The 
middle- left  window  corresponds  to  the  field  evolution,  |o(z,«)|,  over  the  passes  n=500 

passes.  Note  that  the  optical  field  shifts  to  the  left  as  the  centroid  of  the  light  pulse  travels 
slower  than  the  speed  of  light  c.  The  final  optical  pulse  shape  |a(z)|  is  shown  in  the  upper 
left  window.  In  the  specific  simulation  output,  the  optical  field  and  the  power  decay 
accordingly  after  n  ~  150  passes. 
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7  =  5  o.=l  j=o 
2  =10  N  =  36  5^  =0.0001 


Figure  12.  Short  Pulse  Evolution  at  Exact  Synchronism. 

The  middle  lower  window  represents  the  weak  field  spectrum  G(v0)and  is  plotted  for 
reference.  The  window  above  that  shows  the  evolution  of  the  optical  power  spectrum 
P(y,n)  over  n  =  500  passes.  The  final  power  spectrum  P(y )  is  in  the  top- middle 
window.  The  pointed  tick- mark  at  the  top  indicates  the  initial  optical  wavelength  derived 
from  the  resonance  condition.  The  thicker  mark  indicates  the  center  of  the  final  spectrum. 

The  lower-right  window  shows  the  dimensionless  optical  power  P(n)  at  the  end 
of  each  pass.  The  middle  right  window  shows  the  evolution  of  each  of  the  electron  energy 
spectrum  /  (v  ,n )  over  n  =  500  passes.  The  upper- right  window  shows  the  final  electron 
spectrum  after  n  =  500  passes.  Again,  the  pointed  tick- mark  shows  the  initial  phase 
velocity  of  the  electron  beam  at  resonance.  The  thicker  tick- mark  indicates  the  position  of 
the  final  average  electron  phase  velocity  (Av ) . 
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In  exact  synchronism,  the  steady  state  power  is  reduced  to  zero.  The  sbwer  speed 
of  the  light  pulse  centroid  is  overcome  by  reducing  the  path  inside  the  resonator  by  AS  . 
To  do  so,  piezoelectric  crystals  move  the  mirrors  slightly  and  d  =  —2AS  /  NX  .  Figure  13 
shows,  an  FEL  has  the  same  design  characteristics  as  the  FEL  in  Figure  12,  but  with 
d  =  0.06  instead  of  d  =  0  . 


j  = 5  °,=1  d  =  0.06 

g=10  N  = 36  =0.0001 


Figure  13.  Short  Pulse  Evolution  at  Relative  Large  Desynchronism. 


The  introduction  of  a  small  amount  of  desynchronism  d  -  0.06  allows  the  FEL  to 
produce  a  steady  state  power  after  n  =  50  passes. 

As  will  be  seen  in  the  chapters  to  follow,  desynchronism  is  important  to  the 
characteristics  of  pulse  evolution.  By  increasing  desynchronism  d,  less  power  is  obtained 
but  the  FEL  becomes  more  stable.  The  final  optical  pulse  is  longer  than  the  electron  pulse 
and  as  a  result,  the  power  spectrum  is  narrow.  Increasing  the  desynchronism  too  much 
will  cause  the  electro- optical  pulses  to  not  overlap  sufficiently  and  consequently,  the  FEL 
will  not  operate. 
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IV.  SIMULATIONS  OF  THE  100KW  TJNAF  FEL 


A.  INTRODUCTION 

The  Thomas  Jefferson  National  Accelerator  Facility  (TJNAF)  FEL  has 
demonstrated  2  kW  output  power  of  tunable  infrared  radiation  [Ref.  29].  This  is  the 
highest  average  power  level  from  an  FEL  in  the  world  and  was  achieved  with  a  design 
leading  towards  the  further  development  of  get  higher  power.  A  schematic  representation 
of  the  system  is  shown  in  Figure  14. 


In  order  to  operate  at  higher  output  power  levels,  the  TJNAF  FEL  utilizes  electron 
beam  recirculation.  Energy  recirculation  or  “energy  recovery”  is  the  process  by  which  the 
same  electron  beam  is  returned  to  the  accelerator.  After  the  interaction  inside  the  wiggler, 
the  recirculated  electrons  arrive  back  at  the  linear  accelerator  “Linac”  180  degrees  out  of 
phase  so  that  their  energy  is  converted  back  into  RF  power  used  to  accelerate  subsequent 
electrons. 

The  energy  spread  that  a  specific  system  can  recirculate  is  limited.  The  limitations 
come  from  the  induced  RF  phase  shift  upon  arrival  at  the  wiggler,  which  is  equal  to 
optical  desynchronism.  The  change  in  desynchronism  alters  the  FEL  gain  and  the  laser 
output  power.  Consequently,  these  optical  power  changes  result  in  a  change  in  the  energy 
of  the  recirculating  electron  beam,  which  can  potentially  lead  to  additional  electron  beam 
loss  on  apertures  and  a  phase  shift  of  the  decelerated  beam.  The  power  variations,  the 
phase  shifts  and  the  beam  loss  can  change  the  beam- induced  voltage  in  the  accelarator 
cavities.  The  process  can  be  refered  to  as  a  “beam  loading  instability”  [Ref.  31].  If  the  rf 
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control  system  does  not  possess  sufficient  gain  and  bandwidth,  the  whole  recirculating 
FEL  becomes  unstable. 

Another  constraint  on  recirculation  occurs  in  bending  magnets  of  static  ble 
magnetic  fields  are  used  for  the  180  degrees  turns.  The  momentum  of  the  electrons 
spread  resulting  from  the  FEL  interaction  leads  to  a  range  of  bending  radii  that  can  drive 
some  of  the  electron  beam  into  the  wall  of  the  vacuum  pipes  and  cause  the  FEL  to  cease 
lasing.  Overall  ,  there  are  many  benefits  to  recirculation,  such  as  the  enhancement  to  the 
wall-plug  efficiency  and  reduction  of  harmful  radiation  in  the  beam  dump,  so 
recirculation  is  desirable. 

In  the  course  of  developing  megawatt  power  levels,  which  is  required  for  a 
weapon,  it  is  advisable  to  proceed  step  by  step.  Studies  at  TJNAF  have  indicated  that  the 
system  can  be  upgraded  to  operate  at  an  increased  power  of  100  kW.  This  upgrade  has 
been  proposed  with  a  projected  demonstration  scheduled  for  2005  [Ref.  11]. 

To  support  this  goal,  simulations  were  conducted  to  determine  gain,  steady- state 
extraction  efficiency  and  the  induced  energy  spread.  The  input  data  for  simulations  were 
design  parameters  chosen  by  the  TJNAF  scientists  [Ref.  32].  To  achieve  100  kW  average 
power,  the  kinetic  energy  of  the  electron  beam  increased  to  Ee  =  2 10  Me V  and  the  pulse 
repetition  rate  to  £2  =  750  MHz.  The  average  electron  beam  power  is  given  by 

Pe=EJ£ea/c  (4.1) 

where  7  =  270 A  is  the  peak  current  in  the  electron  micropulse,  i e  =0.1  mm  is  the 
electron  micropulse  length  and  c  is  the  speed  of  light.  The  average  electron  beam  power 
then  calculated  as  Pe  =  14  MW.  Therefore  an  output  power  of  100  kW  requires  an 
extraction  efficiency  of  approximately  r)  =  0.7  %  . 

The  undulator  period  is  A,()  =  8  cm  over  N  =  36  periods  with  an  rms  undulator 
parameter  of  K  =  1.7.  Table  1  summarizes  the  experimental  parameters. 
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Electron  energy 

Ee  =  210  MeV 

Peak  electron  pulse  current 

/  =  270  Amp 

Electron  pulse  length 

£e  =  0.1  mm 

Number  of  undulator  periods 

A  =  36 

Undulator  wavelength 

Xa  =  8  cm 

Undulator  length 

L  =  288  cm 

Resonator  length 

S  =  36  m 

Optical  wavelength 

X  =  1  mm 

Step-Taper 

A  =  0,±7t,±27t 

Table  3.  TJNAF  lOOkW  FEL  Parameters. 


In  the  simulations  presented  here,  different  undulator  designs  were  investigated  to 
access  their  possible  advantageous  effects  on  the  FEL  interactions.  The  conventional 
undulator  possesses  a  periodic  magnetic  field  and  wavelength,  but  a  linearly  tapered 
undulator  gradually  changes  the  undulator  parameter  K  as  shown  in  Figure  15  by 
modifying  the  gap  between  the  undulator  magnets  along  the  undulator  length.  The  field 
can  either  increase  or  decrease  along  the  undulator  length. 


LINEAR  TAPER,  8 


8  <  0 


8  >  0 


Magnetic 

Field 


z 

Figure  15. 


Schematic  for  a  Linear  Taper  Wiggler. 
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A  step- tapered  undulator  abruptly  changes  the  value  of  the  field  and  K  halfway 
through  the  undulator  as  Figure  16  shows.  The  stepped  field  can  decrease  or  increase  as 
shown. 

STEP-TAPER,  A 


Magnetic 

Field 


A  >  0 


A  <  0 


ji 

▲ 

titrn 

A 

▲ 

T 

T 

▼ 

T 

▼ 

Figure  16. 


Schematic  for  a  Step-Taper  Wiggler. 


The  resulting  FEL  interaction  in  a  tapered  undulator  [Ref.  25],  [Ref.  33]- [Ref.  37] 
is  described  by  the  modified  pendulum  equation 


V  =f  =  5  +0 


f  n 

T - 

v 


A  +  |a|cos(£  +cp) 


(4.2) 


f  0  for  z  <  0 

where  0  (z)  =  i  is  the  step  function, 

1  for  z  >  0 


C,  =  (k+k0)z-(Ot  is  the  electron  phase,  k  =  2k  / X  is  the  optical  wavenumber, 
k0  =  2k  l\ ,  is  the  undulator  wavenumber,  co  -kc  is  the  optical  frequency,  v  is  the 
electron  phase  velocity,  cp  is  the  optical  phase, 


I  I  ,(P 

a  =  \a\e 


(4.3) 


a 


is  the  complex  dimensionless  optical  field,  |a|  =  2%  NeKL\E\h{  "me2 ,  and, 


6  =  — 4tc  NK~ 


(8  K) 


y 


(1 +K) 


(4.4) 
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where  8  is  the  phase  acceleration  caused  by  the  linear  tapering  of  the  undulator,  b  K  /  K 
is  the  fraction  change  in  the  undulator  parameter  along  the  whole  undulator  length, 


A  =  -4n  NK2 


A K  ' 


(1  +  /T) 


is  the  step- taper  along  the  undulator  at 


T  = 


Z 

NK 


i_ 

2’ 


(4.5) 


(4.6) 


and  A K  IK  is  the  step  change  in  the  undulator  parameter  at  %  =1/2. 

B.  WEAK  FIELD  GAIN 

The  100  kW  TJNAF  FEL  is  described  by  the  dimensionless  current  j  =  5,  the 
electron  pulse  length 


(4.7) 


and  the  resonator  cavity  quality  factor  of  Q  =  4.2  .  To  study  the  FEL  gain,  a  weak  optical 
field  is  used  so  that  the  dimensionless  field  is  \a\  <  n  The  desynchronism  d  is  given  by 


d  = 


AS 

NX 


(4.8) 


and  measures  the  shortening  of  the  resonator  cavity  length  by  AS  compared  to  the 
slippage  distance  NX  .  The  desynchronism  between  the  optical  pulse  and  the  electron 
pulse  was  varied  from  d  =  0  to  d  =0.3.  The  FEL  gain  results  from  many  simulations  at 
different  values  of  desynchronism  d  with  a  step- tapered  undulator  are  plotted  in  Figure 
17. 
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Figure  17.  Weak  Field  Gain  G  versus  Desynchronism  d. 


For  the  conventional  non-tapered  case,  (A  =  0) ,  it  can  be  seen  that  the  highest  gain 
is  60%.  At  small  and  large  values  of  desynchronism,  the  gain  decreases  for  all  the 
undulators  away  from  peak  values  around  d  =0.1.  The  gain  for  small  step- tapers  of 
A  =  ±7l  is  only  slightly  reduced  from  the  conventional  case.  Larger  step-tapers  A  =  ±271 
cause  a  significant  reduction  in  gain  down  to  around  40%,  peaking  at  the  lower  values  of 

d  -0.05.  Also,  the  FEL  still  works  above  the  threshold  losses  (Q  =  4.2)  for  larger 
values  of  d  <  0.2 . 

Figure  18  shows  the  gain  spectrum  for  an  FEL  with  an  inverse  taper  rate  of 
A  =  -27t  over  N  =  36  periods  with  a  current  density  of  /  =  6 .  These  parameters  are 

descriptive  of  the  proposed  TJNAF  100  kW  FEL. 
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ln(lfG) 


Figure  18.  Gain  vs.  v0  for  Inverse  Step  Taper  A  =  -271. 

For  a  negative  step  taper,  the  gain  decreases  in  strong  fields  but  still  remains 
significant.  Also,  increasing  the  optical  field  ao  results  in  a  small  increase  in  the  phase 
velocity  for  peak  gain. 

In  the  electron  pendulum  phase-space,  determined  by  Equation  (4.2),  the  negative 
step  taper  A  =  -2k  FEL  shows  1%  efficiency  with  3.8%  induced  energy  spread.  Figure 
19  shows  a  simulation  in  a  strong  optical  field  where  electron  bunching  is  evident  in 
phase  space  evolution. 

****  pel  Phase  Space  Evolution  **** 
j  =  5  a0=  15  v0=7.4  N  =  36 

Cg=1  A=-2ji 


U j  Sain  0.25 


Figure  19.  Phase  Space  for  Step-Taper  A  =  -2  K  Undulator. 
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As  can  be  seen  in  the  figure  the  electrons  start  at  initial  phase  velocity  v0  =  7.4 
with  a  random  spread  in  phase  velocities  of  width  ov  =  1 .  The  light  colored  dots  show  the 
phase  space  position  at  the  beginning  of  the  undulator,  while  the  darker  colored  dots 
show  the  position  at  the  end.  At  phase  position  C,  ~  n ,  the  electrons  are  bunched,  amplify 
the  optical  field  and  lead  in  a  gain  G  -23%.  The  optical  phase  at  this  position  is  small 
<f>  -0.02  and  further  reduced  to  (f)  -0  at  the  end  of  the  undulator. 

C.  STEADY  STATE  EXTRACTION  EFFICIENCY 

For  a  conventional  undulator  with  no  taper,  the  desynchronism  was  varied  from 
d  =0.005  up  to  0.3  to  study  the  pulse  slippage  effects.  The  best  results  for 
desynchronism  were  d  =0.04,  with  an  efficiency  of  r)  =0.08  above  the  requirement  for 

lOOkW.  At  that  value  of  desynchronism,  the  induced  energy  spread  was 
Ay  /y  =  Av  I  An  N  =  3%,  which  is  well  below  the  15%  limit  for  recirculation.  Figure  20 

shows  that  slight  increases  in  desynchronism  to  d  >  0.06  makes  the  efficiency  marginal 
(r|  <  0.7  %)  for  reaching  100  kW. 


Figure  20.  Efficiency  rj  Versus  Desynchronism  d. 
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Additionally  for  a  step- taper  undulator,  the  results  of  longitudinal  multimode 
simulations  are  shown  in  the  same  figure.  Again,  the  efficiency  is  plotted  versus 
desynchronism  d  for  different  values  of  A  .  For  the  same  value  of  d  =0.04,  a  negative 
step-taper  A  =  -n  once  again  demonstrated  the  best  results,  and  the  efficiency  increased 
to  -0.9%. 

Figure  21  presents  a  high  power  example  for  the  A  =  -7t  case  with  a 
desynchronism  value  of  d  =  0.04  . 

*******  FEL  Puiss  Evolution  ******* 

7=5  <7,  =3  d  =  0.04 

Q  =  4.2  A  =  — 7t  A  =  36 


-6  z  6  -50  V  50  -50  V  50 


-6  z  6  -50  V  50  0  n  500 


Figure  21 .  Power  Simulation  for  Inverse  Step-taper  Undulator. 

The  efficiency  for  this  case  is  r|  =  0.86%  ,  which  corresponds  to  an  output  power  of  116 
kW.  The  induced  electron  energy  spread  is  only  Ay  /y  =  2.4%  ,  which  is  within  the  15% 

specification  for  recirculation.  The  power  reaches  a  steady  state  in  less  than  n  =  50 
passes  and  remains  steady  throughout  the  whole  experiment.  The  good  behavior  of  the 
power  evolution  P(n)  is  indicated  at  the  bottom  right  window  (steady  state  power  after 
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500  passes).  Larger  values  of  tapering  A  <-271 ,  and  desynchronism  d  >0.06  did  not 
work. 

With  both  the  linear  and  step  taper,  a  slight  increase  in  efficiency  over  the 
conventional  undulator  was  observed  for  the  small  negative  taper. 

The  induced  energy  spread  did  not  change  much  with  any  undulator  design  as  is 
shown  in  Figure  22. 


desynchronism  d 


Figure  22.  Energy  Spread  Ay/y  versus  Desynchronism  d  for  Step-Taper  with  Q=4.2. 
Lower  Power  Reduces  Energy  Spread  for  all  Undulator  Designs. 

The  maximum  values  of  the  induced  electron  energy  spread  were  found  in  the  region 
with  maximum  efficiency.  Taper  rates  A  =  0,  -7t  result  in  energy  spreads  of 
Ay/y  =3.5%  and  3.7%  respectively.  These  values  were  the  largest  found  for  all 
undulator  designs  and  were  far  less  than  the  maximum  allowable  limit  for  safe  beam 
recirculation  around  15%.  For  this  reason,  the  quality  factor  Q  was  increased  from  4.2  to 
10.  The  new  resonator  cavity  will  only  have  10%  losses  per  pass.  Also,  larger  efficiency 
and  induced  electron  energy  spreads  were  examined. 
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D.  HIGHER  Q  EXPERIMENTS 

The  purpose  of  these  simulations  was  to  explore  regions  beyond  the  design 
parameters  and  find  the  optimum  undulator  design.  The  quality  factor  Q  was  increased 
from  4.2  to  10  while  the  rest  of  the  input  data  remained  similar.  Such  experiments  have 
been  conducted  at  TJNAF  and  are  described  in  the  following  chapter. 

1.  Step-Taper  Wiggler 

One  of  the  proposed  designs  was  to  make  use  of  the  step-tapered  wiggler.  For  this 
design  the  extraction  efficiency  and  the  induced  electron  energy  spread  were  explored. 
a.  Extraction  Efficiency 

In  the  simulations,  the  step-taper  rates  |A|  were  applied  from  0  — >  471 .  In 
Figure  23  the  curves  show  the  result  of  simulations  with  a  quality  factor  of  Q  -10  . 


Figure  23.  Efficiency  rj  versus  Desynchronism  d  for  Step-Taper,  and  Higher  Q=10. 
Power  Above  1 00  kW  for  a  Larger  Range  of  T apers  A  =  0,±7t,±27t. 


In  this  case,  the  extraction  efficiency  for  some  values  of  A  and  d  exceeds 
1%.  In  this  case  of  A  =  -7T  ,  the  efficiency  increased  from  0.9%  to  1.4%.  For  all  values  of 
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A  =  ±471  and  for  a  desynchronism  of  d  >  1.2,  the  FEL  failed  to  exceed  the  efficiency 
limit.  Compared  to  the  case  of  Q  =  4.2 ,  there  is  a  larger  range  of  desynchronism  in  which 
the  FEL  operates  above  the  extraction  efficiency  requirement  of  r|  >0.7%  .  Moreover, 
other  taper  rates,  such  as  A  =  +7t ,  ±271  ,  were  able  to  produce  output  power  of  exceeding 
100  kW.  This  result  is  significant,  because  it  was  not  confined  to  a  small  region  of 
parameters.  See  Figure  7  for  Q  =  4.2.  More  acceptable  taper  rates  of  A  and  a  larger  d 

lead  to  more  stable,  and  therefore,  more  desirable  systems.  The  simulation  output  for  the 
optimum  inverse  taper  of  A  =  -71  at  d  =  0.03  is  presented  in  Figure  24. 


******* 

FEL  Pulse  Evolution  ******* 
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Figure  24.  Power  Simulation  for  A  =  -;c,<2=10. 

In  this  figure,  the  steady  state  power  is  evident  and  the  induced  energy  spread  is  limited. 
The  efficiency  q  =  13.6%  results  in  an  output  power  of  190  kW  which  is  almost  double 
that  of  the  objective.  Also,  for  this  power,  the  induced  energy  spread  Ay  /y  =  4.3%  is 
again  within  a  safe  limit. 
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b.  Induced  Energy  Spread 

The  new  design  in  the  cavity  Q  - 1 0  increased  the  induced  electron 

energy  spread  significantly.  Nevertheless,  as  Figure  25  shows,  it  is  still  below  the 
maximum  acceptable  limit  of  15%  and  our  design  can  recirculate  the  electron  beam. 


Figure  25.  Energy  Spread  Ay/y  for  Step-Taper,  <2=10. 

2.  Linear  Taper  Wiggler 

A  linearly  tapered  wiggler  gradually  changes  the  undulator  parameter  K  by 
modifying  the  gap  between  the  undulator  magnets.  A  schematic  representation  was  given 
in  Figure  15. 

a.  Extraction  Efficiency 

In  simulations,  the  linear  taper  rate  5  applied  was  varied  over  the  values 
5  =  0, ±7i , ±2% , ±47t ,  and  corresponds  to  AK/ K  =±0.9%, ±1.9%, ±3.8%  .  Recall  from 

equation  (4.4)  that  8  =  -4jc  NK 2  (AK  /  K)(l  +  K2) . 

The  operating  rate  of  desynchronism  d  between  the  optical  and  electron 
pulses  was  varied  from  0  —>0.3.  Figure  26  shows  the  extraction  efficiency  as  a  function 
of  desynchronism  and  the  taper  phase  acceleration  8  . 
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Figure  26.  Efficiency  rj  versus  Desynchronism  d  for  Linear  Taper  with  Q=10.  Power 

Above  100  kW  for  a  Larger  Range  of  Tapers. 

Once  again,  it  is  found  that  the  negative  taper  results  in  maximum 
efficiency.  Lor  the  linear  taper  rate  6  =  -271 ,  there  was  a  maximum  efficiency  of 
r|  -1.4%  at  cl  =  0.01 .  Lor  the  untapered  undulator  6  =  0,  there  was  a  maximum 

efficiency  of  r|  =  1 .2%  at  the  same  value  of  d.  Other  undulators  with  negative  taper  rates 
of  5  --K  ,  —4k  reach  r|  =1.3%,  which  is  better  than  the  zero  tapered  undulator.  Note 
that  1.0%  efficiency  is  easily  achieved  from  almost  all  taper  rate  values  of  8  with  the 
exception  of  positive  tapers  5  =  2tc  ,  4zr  .  Lor  these  taper  values,  r\  =0.9%  and  r\  =0.8% 

occurred  respectively.  These  values  of  extraction  efficiency  are  the  smallest,  but  still 
result  in  an  output  average  power  of  more  than  100  kW. 
b.  Induced  Energy  Spread 

A  smaller  exhaust  energy  spread  allows  the  intense  electron  beam  to 
recirculate  without  damaging  components.  Ligure  27  shows  the  final  induced  energy 
spread  produced  by  seven  different  designs. 


52 


desynchronism  d 


Figure  27.  Energy  Spread  Ay  /y  for  Linear  Taper,  <2=10. 

The  inverse  taper  undulator  for  8  =  -271  has  a  peak  energy  spread  of 
Ay  /y  =  6.3%  occurring  at  d  =0.01.  All  the  values  of  all  taper  rates,  including  the  latter 
value,  produce  an  energy  spread  below  the  maximum  tolerable  limit.  An  increase  in 
desynchronism  in  order  to  reduce  the  electron  beam  energy  spread  is  not  necessary.  Note 
that  a  positive  linear  taper  produces  the  smallest  energy  spread  compared  to  all  other 
designs.  In  small  values  of  desynchronism,  all  the  taper  rate  values  are  almost 
comparable.  However,  after  d  =  0.04 ,  the  values  that  occur  from  the  positive  taper  are 
significantly  smaller. 

E.  CONCLUSIONS 

The  best  results  of  100  kW  simulations  for  various  designs  and  a  quality  factor  of 
Q  =  4.2  are  summarized  below  in  Table  4. 
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2  =  4.2 ,  STEP  TAPER  DESIGN 

Step-Taper 

Desynchronism 

Efficiency 

Induced  energy 

Output  Power 

A 

d 

r|(%) 

spread  Ay/y  (%) 

P  (kW) 

0 

0.04 

0.82 

3.5 

115 

K 

0.04 

0.55 

3.2 

77 

2n 

0.05 

0.25 

2.4 

35 

-K 

0.04 

0.86 

-3.5 

120 

-2% 

0.05 

0.64 

3.2 

90 

Table  4.  Results  of  TJNAF  100  kW  Step-Taper  Undulator,  for  Q  -  4.2 . 

Step-taper  wigglers  were  explored  and  compared  with  conventional  wigglers.  For 
the  given  TJNAF  lOOkW  FEL  parameters,  the  inverse  taper  undulator  A  =  -n  achieved 
the  highest  final  power  at  a  desynchronism  of  d  =  0.04  .  The  output  power  produced  was 
P  -  1 20  kW  which  is  above  the  100  kW  objective.  The  conventional  untapered 
undulator,  (5  =0),  produced  a  final  power  of  P  =115  kW,  which  again  was  above  the 
objective  at  the  same  desynchronism  value.  In  both  cases,  the  induced  electron  energy 
spread  Ay  /y  -  3.5%  was  smaller  than  the  maximum  allowable  limit  of  15%  and  did  not 
limit  the  operational  range  of  the  FEL. 

The  small  values  of  Ay  /y  allowed  for  the  possibility  to  explore  a  design  with 
10%  losses  (Q  =10)  in  an  optical  cavity.  The  high  Q  =10  factor  was  chosen  in  such  a 
way  as  to  increase  the  efficiency  of  the  system,  while  maintaining  the  induced  energy 
spread  well  below  the  15%  engineering  limit.  The  results  related  with  Q,  this  design  and 
the  step-taper  design  are  summarized  in  Table  5. 
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2=10,  STEP  TAPER  DESIGN 


Step-Taper 

Desynchronism 

Efficiency 

Induced  energy 

Output  Power 

A 

d 

r|(%)  " 

spread  Ay/y(%) 

P  (kW) 

0 

0.01 

1.19 

5.3 

167 

K 

0.01 

0.94 

5.1 

132 

2k 

0.01 

0.78 

5.6 

109 

4k 

0.02 

0.12 

1.5 

17 

-n 

0.01 

1.36 

5.1 

190 

-2k 

0.01 

1.25 

5.3 

175 

-4k 

0.03 

0.45 

3.5 

63 

Table  5.  Results  for  Step-Taper,  Q  =  1 0  Design. 

The  best  result  for  all  the  above  cases  was  a  taper  rate  of  A  =  -n  and  occurred  at 
desynchronism  d  =  0.01  with  an  output  power  of  P  =  190  kW.  Moreover,  all  the  designs, 

except  A  =  ±471  ,  operated  above  the  threshold  and  are  capable  of  producing  an  output 
power  of  at  least  100  kW.  The  linear  taper  results  are  summarized  in  Table  6. 


2=10,  LINEAR  TAPER  DESIGN 


Linear  Taper 

Desynchronism 

Efficiency 

Induced  energy 

Output  Power 

8 

d 

Tl(%) 

spread  Ay/y (%) 

P  (kW) 

0 

0.01 

1.19 

5.3 

167 

K 

0.01 

1.1 

5.5 

154 

2k 

0.01 

0.92 

5.1 

129 

4k 

0.01 

0.68 

5.5 

95 

-K 

0.01 

1.31 

5 

183 

—2k 

0.01 

1.43 

5.2 

200 

-4k 

0.01 

1.32 

5.1 

185 

Table  6.  Linear  Taper,  2=10  Design. 
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The  linear  taper  seems  to  be  a  very  promising  design.  All  taper  rates,  except  for  5  =  471  , 
achieved  100  kW.  For  5  =-271  at  a  desynchronism  of  d  =0.01,  the  output  power 
P  =  200  kW  was  double  the  objective.  From  the  table,  the  superiority  of  the  inverse  taper 
is  clear.  This  kind  of  taper  results  in  a  large  efficiency  while  the  induced  electron  spread 
is  comparable  to  other  designs,  and  sufficiently  lower  than  the  15%  limit  for  safe 
recirculation. 

The  common  behavior  for  all  forms  of  taper  rates  is  a  reduction  in  extraction 
efficiency  and  the  energy  spread  for  a  slight  increase  in  desynchronism. 

This  work  was  presented  at  the  23rd  International  FEL  Conference  in  Darmstadt, 
Germany  and  has  been  accepted  for  publication  in  Nuclear  Instruments  and  Methods  in 
Physics  Research. 
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V.  THE  TJNAF  TAPERED  WIGGLER  EXPERIMENT 


A.  INTRODUCTION 

Using  the  Thomas  Jefferson  National  Accelerator  Facility  (TJNAF)  free-electron 
laser,  experiments  have  been  conducted  which  study  its  operation  as  a  function  of  various 
taper  rates  and  desynchronism  values  [Ref.  38].  The  kinetic  energy  of  the  electron  beam 
was  Ee  =  34.5  MeV,  the  optical  wavelength  was  X  =  6  |im,  the  electron  pulse  duration 

was  0.5  ps  (pulse  length  i e  =150pm),  and  the  peak  current  was  I  =50A.  The  purpose 
of  the  simulations  presented  here  is  to  compare  experiment  with  theory,  and,  by  using  a 
wider  range  of  parameters  than  allowed  by  experiment,  to  extract  more  physical  meaning 
from  the  results. 

B.  PARAMETERS  AND  SIMULATIONS  METHODS 

The  simulations  used  the  same  parameters  as  TJNAF,  described  by  dimensionless 
current  density  /  =  1 0  ,  normalized  electron  pulse  length  a  =  NX  =  1  and  undulator 

parameter  K  =0.98(rms).  Table  7  summarizes  the  experimental  parameters  and  Table  8 
the  corresponding  dimensionless  parameters  that  are  used  in  the  simulations.  The 
desynchronism  between  the  optical  pulse  and  the  electron  pulse  was  varied  from  d  =  0 
up  to  the  highest  value  of  d  which  produced  power  corresponding  to  an  efficiency  p  > 
0.01%.  As  an  example,  the  maximum  desynchronism  for  the  conventional  non-tapered 
case  (8=0)  is  d  =  0.4 .  Not  only  were  taper  values  for  which  we  had  experimental 
results  explored,  i.e.  5=0,  10rc ,  -In  ,  but  values  well  beyond  these  were  also  explored. 
By  using  a  wider  range  of  values  of  linear  taper,  it  was  possible  to  plot  the  operating 
width  of  the  desynchronism  curve  versus  linear  taper.  The  efficiency  p  and  the  induced 
energy  spread  Ay  /y  are  determined  as  functions  of  desynchronism  d  and  taper 
acceleration  5  . 
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TJNAF  EXPERIMENTAL  PARAMETERS 


Electron  Energy 

y  me2  =34.5  Mev 

Undulator  Length 

L=  1.1  m 

Undulator  Periods 

N=  41 

Undulator  Wavelength 

Xo  =2.7  cm 

Peak  Undulator  Magnetic  Field 

B  =5.5  kG 

Peak  Current 

7  =  50  A 

Taper  Phase  Acceleration 

(8/7t)  =-7,0,  10 

Cavity  Losses 

Q  =50 

Table  7.  Parameters  used  in  the  Experiment 


DIMENSIONLESS  SIMULATION  PARAMETERS 

Current  density 

7=10 

Electron  Pulse  width 

<u=i 

Undulator  Parameter 

K=  0.98 

De  synchronism 

d=  0  —>0.4 

Taper  Rate 

(8  In )  =  0,  ±3,  ±6,  ±7,  ±8,  ±9,  ±  10,  ±12,  ±  14,  ±16,  ±  18,  ±20, ±22,  ±24 

Table  8.  Parameters  used  in  the  Simulations 
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C.  UNTAPERED  UNDULATOR 

For  an  untapered  undulator,  the  desynchronism  was  varied  from  d  =  0.002  up  to 
0.4  in  order  to  study  pulse  slippage  effects.  It  was  found  that  for  desynchronism 
d  =0.02,  the  efficiency  was  r|=0.04.  However,  at  this  value  the  power  P(n)  is  not 
steady,  as  shown  in  Figure  28  and  the  power  oscillations  are  caused  by  the  trapped- 
particle  instability  resulting  from  strong  optical  fields  [Ref.  25]. 

****  pgL  Pulse  Evolution  **** 
j  =  10  oz=l  d  =  0.02 
Q  =  50  N  =  41  8C,  =  0.0001 


Figure  28.  Trapped  Particle  Instability. 


In  the  upper  left  window  the  final  optical  pulse  shape  after  n  =  500  passes  is 
shown.  There,  short  optical  pulses,  modulated  with  sharp  spikes,  create  two  sidebands  in 
the  optical  power  spectrum  P(v,n)  shown  in  the  upper- middle  plot.  In  the  upper  right 
window,  the  electron  spectrum  f(v,n)  shows  the  large  induced  spread  in  phase  velocities 
Av  =44.  which  is  not  desirable  for  recirculation  of  the  electron  beam.  The 
corresponding  induced  energy  spread 


Ay  _  Av 

y  AnN 


44 
47t  N 


=  9% 


(5.1) 
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Steady-state  power  without  the  trapped-particle  instability  or  limit-cycle  behavior 
[Ref.  39]  is  shown  in  Figure  29  for  larger  desynchronism  of d  >0.15.  Note  the  narrow 
electron  distribution  and  the  absence  of  sidebands  in  the  power.  Also  note  that  the  power 
evolution  P(n)  (bottom  right)  has  no  oscillations. 


****  FEL  Pulse  Evolution  **** 
j  =  10  Cfz  =1  d  =  0.15 
Q  =  50  N  =  41  5^  =  0.0001 


Figure  29.  Steady  State  Power  for  Conventional  Undulator. 


D.  GAIN  SPECTRUM 

The  dimensionless  electron  phase  velocity  v  .  is  the  meeting  point  between 
simulations  and  experiments  and  in  dimensionless  notation  is  given  by 


v  =  2nN 


(i+k2) 

2y2l 


K 


(5.2) 


where  N  =  41  is  the  number  of  undulator  periods  and  y  is  the  relativistic  Lorentz  factor. 
Single  mode  simulations  are  used  to  find  the  optimum  initial  phase  velocity  vo  for 
maximum  gain. 
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Figure  30  presents  a  plot  of  vo  as  a  function  of  the  tapering  rate  5  .  The  relation 
for  peak  gain  vo  =  2.4-8  / 2  is  presented  with  the  red  circles  and  is  valid  in  weak  fields 
and  moderate  taper.  Taper  rates  8  >1071  and  strong  fields  \a0  \  >n  change  this  relation. 


Figure  30.  Plot  of  Phase  Velocity  versus  Tapering  Rate. 

Figure  31  presents  the  peak  gain  as  a  function  of  taper  by  using  the  optimum 
initial  phase  velocity  vo  of  Figure  5.  Is  required  for  all  tapering  rates  -2471  <8  <  24n, 

we  have  gain  G  that  exceeds  the  resonator  losses.  The  optical  cavity  losses,  used  in  the 
experiment,  were  for  Q  =  50  so  that,  gain  G  >  2%  is  required  for  operation.  The 

minimum  gain,  after  introducing  taper  was  G  =  7% ,  and  was  met  for  8  =  2471  . 
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The  G(S)  curve  is  almost  symmetric  in  8  with  peak  value  at  8  =  0  where  the  gain 
is  G=160%.  As  the  magnitude  of  8  increases,  the  peak  gain  decreases  and  for 
8  =  ±2471  reduces  to  G  ~  8%  .  There  is  only  a  1%  difference  in  gain  at  the  ends  of  the 
curve  where  negative  taper  results  in  slightly  more  gain. 

Figure  32  (a)  shows  single  mode  gain  plotted  as  a  function  of  the  initial  optical 
field  ao  and  the  initial  electron  phase  velocity  vo.  For  8=0,  weak  fields  (a0<  7t)  and 

current  j= 10,  the  maximum  gain  is  G=160%  and  occurs  at  v  =2.4.  Figure  32  (b) 

shows  that  the  application  of  modest  positive  taper  rate  (8  >0)  makes  the  peak  gain 

decrease  and  shift  to  the  left  in  va.  Figure  32  (c)  shows  similar  results  for  the  inverse 

taper  8  <0but  shifts  to  a  reduced  peak  gain  to  the  right  in  va.  For  both  positive  and 

inverse  tapers,  it  can  be  seen  that  gain  is  not  symmetric  around  the  resonance  (v0  =0)  and 

for  large  values  of  optical  field  strength  ao,  there  is  a  distortion  in  the  shape  of  the  curve. 
Mention  that  positive  taper  has  one  gain  peak  in  strong  fields  while  negative  taper  has 
two,  nearly  equal,  gain  peaks  in  strong  fields. 
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In (1+G) 


Figure  32. 


In (1+G) 


j=10  N=41 


0.44 


In  (1+G) 

A 

0.43 

j=10  N-41 

5=-8x  a=0  » 

mth 


(b) 


(c) 

G(v0,a0)  for  Zero  taper  (a),  8  =8ti  (b),  8  =  -871  (c). 


63 


E.  STEADY  STATE  POWER 

In  our  multi- mode  simulations,  the  following  values  of  linear  taper  were  used: 

5  =  ±3tc  ,±6rc  ,±771+871  ,±9Jt,±107t, ±1271, ±147t,±167t, ±187t,±207t, ±2271, ±247t. 

Figure  33  summarizes  the  results  of  multimode  simulations  and  plots  the 
efficiency  rj  for  positive  taper.  The  untapered  undulator  reaches  its  maximum 
efficiencyr|  =  4%  at  d  =0.002  .  As  the  taper  rate  increases  the  efficiency  decreases.  For 
8  =9jt,  the  efficiency  is  50%  smaller  and  peaks  at  d  =0.006.  For  5  =14tc,  the 
efficiency  diminishes  significantly  to  r|  =0.8%  at  d  =0.006  .  After  5  =  14tc  ,  the  peaks  in 
the  efficiency  curves  stay  at  lower  values  of  desynchronism.  Sidebands  are  formed  due  to 
the  trapped-particle  instability  and  the  efficiency  is  enhanced  up  to  2%  for  5  =  1671  at 
d  =0.002.  Finally  for  larger  taper  rates  5  =  247t,the  efficiency  drops  to  only  1%  at 

d  =0.002. 


Figure  33.  Efficiency  rj(d,8)  for  Positive  Taper. 


64 


In  Figure  34  the  efficiency  is  plotted  as  a  function  of  desynchronism  d  for 
negative  taper  rates.  As  the  taper  rate  increases  from  5  =  0  — >  —  3xc ,  the  efficiency 
decreases  from4%— >  2.6%  .  The  peak  of  the  efficiency  curve  continues  to  stay  at  small 
values  of  desynchronism. 


Figure  34.  Efficiency  r|(d,8)  for  Negative  Taper. 


For  modest  taper  rates  5  =  -6ji  ,  -Ik ,  -871 ,  the  efficiency  increases  from 
T|  =  2.9%  — >  3.3%.  With  taper  rate  of  5  =  -9k  ,  maximum  efficiencyr|  =  4.2%  is  reached 
at  d  =  0.002 .  Again  sidebands  are  present  and  responsible  for  this  enhancement  in 
efficiency.  For  5  =  —  1  4tt  — >  -20rt  the  peak  efficiency  is  shifting  slightly  from 
d  =  0.002  — >  0.004  and  efficiency  diminishes  from  r|  =3.3%— >1.7%  .  Further  increase 

of  the  taper  rate  5  — >-247t,  results  in  a  decrease  in  efficiency  to  1%  at  slightly  larger 
d  =0.006. 

Figure  35  shows  the  width  of  desynchronism  Ad  where  the  FEL  can  operate. 
Values  of  Ad  against  the  taper  rate  8  were  plotted. 
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Figure  35.  Effective  Desynchronism  Width. 


Note  that  Ac/(S)  has  the  same  shape  as  the  gain  G(S)  curve  in  Figure  31,  with  a  maximum 
at  8  =  0  and  decreasing  rapidly  as  |S  |  increases.  In  the  region  where  |S  |  ~  ±871  the  rate  of 
descent  becomes  much  less  rapid,  indicating  that  something  special  happens  there. 

The  change  in  slope  near  8  «  ±871  is  explained  by  looking  at  the  gain  curve  (gain 
versus  phase  velocityv0)  as  8  changes.  Figure  36  shows  that  there  are  two  gain  peaks  of 

equal  height  8  =  -871  .  There  are  also  two  gain  peaks  when  8  =  +87t ,  but  on  the  opposite 

side  of  resonance  (v  =  0)  where  (v  <  0).  These  peaks  act  together,  in  either  the  8  =  ±871 

cases,  to  further  amplify  the  optical  field.  The  transition  in  plateau  at  8  =  ±871  results 
from  the  FEE  switching  for  one  gain  peak  to  a  second. 
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0.2431 


0.0 


-0.2431 

0.16 


0.0 


-0.16 


Figure  36.  Comparable  Peaks  in  Gain  Spectrum,  8  =-8rc. 

F.  ENERGY  SPREAD 

The  energy  spread  resulting  from  the  FEL  interaction  can  potentially  limit  the 
performance  of  a  high  power  FEL  using  energy  recovery.  Energy  recovery  is  desirable 
because  it  can  significantly  increase  the  FEL  efficiency  and  reduce  the  size  of  beam 
dump.  In  the  experiment  we  describe,  electron  beam  recirculation  is  considered  feasible 
only  if  the  induced  energy  spread,  Ay/yis  less  than  6%.  If  it  is  larger,  then  bending 

magnets  are  not  able  to  feed  back  the  beam  to  the  accelerator  and  the  FEL  stops 
operating. 

In  order  to  take  into  account  that  limiting  factor,  we  measure  the  width  Av  e  from 
the  electron  distribution  f(y,n)  in  the  FEL  pulse  evolution  simulations.  Accordingly  we 
calculate  the  energy  spread  using  Ay  /y  =  Av^  /  4tc  N  ,  where  N  =  41 .  In  Figure  37,  the 
energy  spread  induced  by  the  FEL  interaction  is  plotted  as  a  function  of  desynchronism  d 
and  positive  taper  8 . 
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The  induced  electron  beam  energy  spread  for  the  conventional  undulator  (5=0) 

has  a  peak  energy  spread  of  14%  occurred  at  d  =  0.002 .  The  strong  fields  cause  the 
trapped- particle  instability  and  consequently  the  energy  spread  is  moved  outside  the  6% 
design  tolerance.  The  induced  energy  spread  decreases  with  increasing  the 
desynchronism,  and  falls  below  the  desired  level  at  a  desynchronism  \aluc  of  d  =  0.04 , 
and  remains  below  the  6%  for  all  larger  values  of  d. 

For  the  positive  tapered  undulator,  the  induced  energy  spread  is  slightly 
decreased.  For  some  values  of  desynchronism  d  <  0.05  and  taper  5  <  10jt  we  are  still 
outside  of  the  limit  that  imposed  for  safe  recirculation.  Unfortunately  this  is  the  region  in 
which  we  have  peak  efficiency.  A  further  increase  in  the  taper  rate  5  >1071  results  a 
smaller  value  for  d  for  which  the  energy  spread  is  acceptable.  The  minimum 
desynchronism  that  reduces  the  energy  spread  down  to  the  desired  6%  is  given  at  the 
second  column  of  the  Table  9. 
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Figure  38  shows  the  final  induced  electron  beam  energy  spread  produced  by 
inversely  tapered  undulators.  The  energy  spread  Ay/y,  as  a  function  of  desynchronism  d 
and  tapering  rate  8,  is  plotted. 


Figure  38.  Electron  Energy  Spread  for  Negative  Taper. 


The  electron  energy  spread  curves  follow  the  same  trend  as  the  power  curves.  For 
high  power,  there  is  a  large  energy  spread.  For  taper  5  =  -971  the  maximum  value  is 
18%  at  d  =0.002.  Recall  that  for  8  =0,  the  maximum  value  is  14%  at  d  =0.002. 
After  d  =0.05  the  energy  spread  drops  to  less  than  6%  for  all  the  tapering  rates  and 
there  is  no  constraint  to  the  FEL  operation.  The  inverse  tapered  undulator  does  not  appear 
to  be  desirable  for  recirculating  the  electron  beam.  This  is  the  trade-off  of  the 
enhancement  in  the  efficiency  with  the  inverse  taper.  The  minimum  desynchronism  for 
which  the  energy  spread  constraint  is  not  met  for  the  inverse  tapered  undulators  is  given 
at  the  fifth  column  of  the  Table  9. 
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Table  9  is  obtained  from  previous  Figures  and  helps  to  summarize  the  values  of 
desynchronism  d  and  taper  8,  where  we  have  acceptable  results  that  meet  the 
requirements. 

•  min  d  for  (Ay  / y  <  6%,  acceptable  energy  spread) 

•  max  d  for  T|  >  0  (FEL  operation) 

In  this  way  we  are  able  to  determine  the  operational  range  of  the  FEL  as  a  function  of  d 
and  8. 


Negative  8 

Min  d 

Max  d 

-371 

0.05 

0.36 

-671 

0.05 

0.22 

-771 

0.05 

0.15 

-871 

0.05 

0.1 

-971 

0.05 

0.08 

-IOtt 

0.05 

0.08 

-1271 

0.05 

0.08 

1471 

0.05 

0.06 

-1671 

0.03 

0.05 

-1871 

0.02 

0.04 

-2071 

0.01 

0.03 

-2271 

0.02 

0.03 

-2471 

0.02 

0.03 

Positive  8 

Min  d 

Max  d 

0 

0.004 

0.4 

371 

0.05 

0.36 

6ti 

0.05 

0.22 

771 

0.05 

0.16 

871 

0.05 

0.1 

971 

0.05 

0.09 

1  0k 

0.03 

0.09 

1271 

0.01 

0.08 

1471 

0.02 

0.06 

1671 

0.02 

0.05 

1871 

0.01 

0.04 

2071 

0.01 

0.03 

22ti 

0.02 

0.03 

2471 

0.006 

0.008 

Table  9.  Desynchronism  Values  Meeting  the  Energy  Spread  Limitation  and  Efficiency. 
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G.  EXPERIMENTAL  AND  SIMULATION  RESULTS  -COMPARISON 

In  the  experiment,  tapers  were  obtained  by  inserting  precision  shims  at  either  end 
of  the  wiggler.  Dial  gauges  on  either  side  of  the  wiggler  measured  the  position  and  gap  of 
the  wiggler  for  each  taper  value.  Since  the  variation  of  the  wiggler  field  with  gap  size  is 
known,  the  field  taper  can  be  calculated  from  the  gap  taper  [Ref.  40] . 

Three  values  of  taper,  8  =  — 7rc  ,0, 1  Otc  were  examined  with  the  following 
experimental  results  compared  to  simulations. 


Efficiency  T| 


Experiment 

Simulations 

8=-7n 

2.0% 

1.54% 

8=0 

1.5% 

1.0% 

S=10jt 

0.62% 

0.67% 

Energy  Spread  Ay/y 


Experiment 

Simulations 

S=-7  n 

2.3% 

3.1% 

S=10jt 

1.6% 

1.8% 

•  Desynchronism  Width  d 


Experiment 

Simulations 

8=-7n 

0.15 

0.10 

8=0 

0.2 

0.35 

S=10jt 

0.075 

0.07 

Table  10.  Experiment  and  Simulations  Results. 

We  used  the  value  of  efficiency  that  allows  safe  recirculation.  For  inverse  taper  5  =  -7n  , 
we  disregard  the  highest  efficiency  r\  =3.2%  at  <7=0.004  because  it  gives  energy  spread 
Ay  /y  -14%  ,  which  is  almost  twice  the  allowed  limit.  The  same  tactic  was  followed  for 
positive  taper  and  so  we  did  not  accept  T|  =1.9%  at  <7=0.002  since  Ay/y  -15%  .  For  the 
zero  taper  case,  the  efficiency  was  lower  than  the  experiment.  The  peak  efficiency  of 
r)  =  3.9%  ,  at  d  =  0.002  ,  was  accompanied  by  energy  spread  Ay/y  - 14%  .  It  is  possible 

to  reduce  the  energy  spread  by  increasing  desynchronism  <7.  At  d  =  0.02  the  efficiency  is 
r)  =  1.54%,  while  the  energy  spread  is  Ay  /y  =  7.5%  .  At  d  =0.03  the  energy  spread  is 
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Ay  /y  =  5.7%  with  r|=1.32%.  For  all  the  above  cases,  safe  recirculation  near 
Ay  /y  ~  6%  was  the  goal.  The  FEL  interaction  worked  well  above  the  engineering  design 
for  the  specific  experiment. 

For  the  energy  spread  Ay  /y  =  Av  /4tt  N  we  used  Av  from  simulations,  when  the 
desynchronism  cl  was  at  the  upper  limit  and  the  efficiency  was  marginal.  In  that  region 
the  optical  and  the  electron  pulse  fail  to  overlap  and  power  was  reduced  to  zero.  For  the 
inverse  taper  5  =-771  that  happened  at  d  =0.1,  where  the  efficiency  was  r)  =0.47% 

and  for  positive  taper  5  =  10ji  that  condition  was  met  at  d  =  0.05  with  r)  =  0.6%  . 

H.  CONCLUSIONS 

For  small  values  of  taper  and  with  small  values  of  cl,  the  FEL  exhibits  the  trapped- 
particle  instability.  Increasing  either  desynchronism  d  or  taper  rate  8  the  FEL  easily 
moves  to  a  more  stable  regime,  even  though  the  efficiency  is  lower.  This  stability  is 
needed  in  some  applications,  where  the  absence  of  side  bands  is  crucial.  For  example, 
the  use  of  the  FEL  in  human  surgery  requires  a  narrow  energy  spectrum  in  order  to 
achieve  the  exact  desirable  wavelength  ( =  6.45  |im)  [Ref.  41]. 

Negative  taper  produces  better  extraction  efficiency  than  positive  taper  as 
theoretical  studies  by  Saldin,  Schneidmiller  and  Yurkov  have  shown  [Ref.  37]. 
Unfortunately  this  is  not  accompanied  by  reduction  in  the  energy  spread  as  T|  is  directly 
related  to  Ay/y  ,  recall  r\  =<Ay me2  >/ymc2 .  At  5  =  -971 .  we  have  a  0.3%  increase 
from  zero  taper,  which  results  r|  =  4.2% .  In  this  manner,  we  not  only  enhance  the 
efficiency,  but  move  away  from  the  trapped-particle  instability  region. 

When  the  FEL  is  operated  with  higher  Q  =  50  ,  the  gain  per  pass  is  significantly 
more  than  the  loss  per  pass.  One  result  was  that  the  range  of  desynchronism  values  when 
the  FEL  operates  with  a  gain  above  threshold  increased  compared  to  the  lower  Q  case. 
Figure  39  shows  the  desychronism  curve  width  Ad  plotted  against  the  positive  and 
negative  taper  rate  8  .  The  superimposed  G(S)  shows  the  clear  correspondence  between 

increased  gain  G  and  increased  desychronism  curve  width  Ad  . 
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Figure  39.  Gain  G(8)  and  Desynchronism  Width  Ad(S). 

The  maximum  gain  occurs  with  no  tapering,  8=0,  as  would  be  expected.  As  the 

taper  rate  8  increases  or  decreases,  the  gain  decreases  from  160%  to  about  40% 
symmetrically  around  8=0.  At  positive  or  negative  taper  rates  of  about  8  =  ±871  ,  the 
gain  plateaus  decrease  much  more  slowly  with  the  increasing  taper  rate.  While  the  taper 
rate  increases  in  magnitude  from  |S  |  =  8tc  out  to  |S  |  =  2471  ,  the  gain  decreases  from  about 
40%  down  to  about  10%. 

The  gain  plateau  at  large  values  of  taper  begins  at  a  value  of  taper  8  that  causes 
the  gain  spectrum  to  change  shape  so  that  there  are  two  peaks  with  nearly  the  same  peak 
gain.  As  the  taper  is  increased  and  the  peak  available  gain  decreases,  the  gain  spectrum 
acquires  multiple  peaks  with  comparable  gain.  As  the  primary  gain  peak  decreases  other 
surrounding  gain  peaks  increase  in  comparison.  As  the  peak  available  gain  increases,  or 
decreases,  at  any  value  of  tapering,  the  desynchronism  curve  width  A d  increases  or 
decreases  correspondingly.  With  more  gain  available  there  are  more  values  of 
desynchronism  that  are  above  the  threshold  making  the  desynchronism  curve  wider. 

The  simulations  and  the  experimental  data  are  in  good  agreement  with  analytic 
theory.  Theory  indicates  that  as  taper  increases,  the  desynchronism  width  decreases.  The 
experiment  shows  this  trend,  but  not  as  clearly  as  the  simulation  data.  There  is  only  one 
significant  deviation  between  simulation  and  experimental  results.  The  desynchronism 
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width  for  zero  taper  is  A  d  -  0.35,  compared  to  Ad  =  0.2  in  the  experiment.  The 

disagreement  has  led  to  reconsideration  of  the  experimental  procedure  for  this  data  point 
[Ref.  41]. 

The  wiggler  design  can  be  important  in  the  energy  recovery  FELs.  The  specific  IR 
TJNAF  FEL  demo  was  designed  to  operate  in  the  safe  mode  where  Ay  /y  <  6%  ,  even  if 
it  were  capable  of  higher  extraction  efficiency  Current  developments  at  TJNAF  make 
feasible  the  recirculation  of  the  electron  beam  for  all  values  of  Ay  /y  that  we  found  in 
the  simulations  runs. 
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VL  CONCLUSIONS 


The  missile  defense  capabilities  currently  under  development  must  not  only 
review  past  wars  but  more  importantly  plan  for  the  future.  Tactical  missile  threats  still  in 
existence  from  past  wars  and  their  capabilities  are  rapidly  improving  with  the  result  that  it 
will  be  much  more  difficult  to  counter  such  threats  in  the  future.  The  current  development 
of  weapons  is  centered  on  shortsighted  decisions  based  on  current  threats  and  technical 
challenges.  By  the  time  some  weapons  are  deployed  they  are  outdated  and  ineffective. 
Light-based  weapons  offer  the  solution  to  this  problem  since  they  possess  layered 
capability  to  counter  a  potential  threat. 

FEL  defensive  capability  prompted  the  U.S.  Congress  to  approve  a  $15  million 
grant  to  the  Thomas  Jefferson  National  Accelerator  Facility.  At  present,  the  TJNAF  FEL 
is  the  most  powerful  FEL  in  the  world  and  operates  at  more  than  2  kW  average  power  for 
many  hours  at  a  time.  One  transitional  step  for  the  1  MW  output  power,  required  for  a 
weapon,  is  the  proposed  100  kW  upgrade.  Multimode  simulations  used  at  NPS  show  that 
this  development  is  feasible. 

For  the  TJNAF  lOOkW  FEL  parameters,  the  inverse  taper  undulator  of  A  =  -7t 
achieved  the  highest  final  power  at  a  desynchronism  of  d  =  0.04  for  Q  =  4.2 .  The  output 

power  produced  was  P  -  1 20  kW,  which  is  above  the  100  kW  objective.  The  small 
values  of  the  induced  electron  energy  spread  (-3.5%)  made  it  possible  to  explore  a  design 
with  10%  losses  (Q  =  10)  in  an  optical  cavity.  The  higher  Q  factor  was  chosen  to 
increase  the  efficiency  of  the  system,  while  maintaining  the  induced  energy  spread  well 
below  the  15%  engineering  limit.  The  results  related  with  the  <2=10  design  were 

•  the  inverse  step- taper  undulator  A  =  -n  achieved  the  highest  final  power 
P  =  190  kW  at  a  desynchronism  of  d  =  0.01 . 

•  The  inverse  linear  taper  5  =  -2n  achieved  the  highest  final  power 
P  =  200  kW  at  a  desynchronism  d  =  0.01 . 

The  validity  of  our  results  have  been  compared  to  experiments  conducted  at  the 
TJNAF  FEL  facility.  Various  taper  rates  and  desynchronism  values  were  studied.  The 
purpose  of  the  simulations  was  to  compare  experiments  with  theory,  and,  by  using  a 
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wider  range  of  parameters  than  allowed  in  the  experiment,  to  extract  more  physical 
meaning  from  the  results.  The  simulations  and  the  experimental  data  were  consistent  with 
analytic  theory.  Theory  indicates  that  as  a  taper  increases,  the  desynchronism  width 
decreases.  The  experiment  shows  this  trend,  but  not  as  clearly  as  the  simulation  data. 
There  is  only  one  significant  discrepancy  between  simulation  and  experimental  results, 
that  has  resulted  in  ongoing  discussion  about  the  validity  of  the  experimental  observation 
in  this  case. 

Thus  it  is  possible  to  cone  lude  that  FELs  have  the  potential  to  be  developed  into  a 
compact,  reliable  and  efficient  weapon.  This  development  leads  not  only  to  a  new 
military  era  but  opens  the  door  for  commercial  and  scientific  applications.  The  extremely 
high  industrial  and  technolo  gical  momentum  that  exists  makes  FEL  development  feasible 
with  high  level  pay-offs. 
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